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ABSTRACT 


A  data  base  is  being  developed  for  calculation  of  quasibinary  and 

quasiternary  phase  diagrams  of  oxide  and  oxynitride  ceramic  systems. 

The  current  base  covers  combinations  of  ^2^1*  MgO>  A^Oj, 

FejO^,  FeO,  CaO,  SiC^,  Si^N^  and  AIN.  Lattice  Stability,  Solution 

and  Compound  Phase  Parameters  were  derived  covering  the  liquid, 

spinel,  periclase,  corundum  alpha  and  beta  quartz,  crystobalite , 

tridymite  and  beta  Si^N^  structures  in  all  of  the  quasibinary 

systems  composed  of  the  data  base  compounds.  Evaluation  of  these 

parameters  permit  calculation  of  the  thermochemical  properties 

of  these  systems  in  good  agreement  with  limited  thermochemical 

data.  In  addition,  isothermal  quasiternary  systems  of  selected 

oxide  and  oxynitride  systems  composed  of  data  base  compounds 

have  also  been  calculated  over  a  wide  range  of  temperature. 

The  free  energy  of  formation  of  beta.sialons  in  the  SijN^A^Oj-AlN 

system  and  the  X  phase  in  the  SiC^-A^O^-Si^N^  system  have  been 

calculated.  Recent  examples  of  the  use  of  data  bases  and  computer 

techniques  in  solving  problems  associated  with:  in-situ  formation 

of  columbium,  nickel  and  cobalt  based  superalloys,  sigma  phase 

formation,  multicomponent  electroptical  semiconductors,  sialon 

ceramics,  molten  salt,  sulfide  and  oxide  slags,  high  temperature 

battery  electrodes,  hardenability ,  synthesis  of  hydrogen  storage 

compounds,  grain  refinement  in  aluminum  alloys,  metatectic  reactions 

in  hafnium  alloys  and  the  growth  of  HEM  silicon  solar  cell  materials. 

The  initial  results  obtained  by  applying  these  techniques  to  CVD 

3 

synthesis  of  doped  germanium-borosilicate  glasses  for  C  I 
applications  are  reported. 
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Section  I.  Introduction  and  Summary 

The  utility  of  computer  based  methods  for  developing 
and  utilizing  coupled  phase  diagram  and  thermochemical  data  on 
metallic  systems  has  been  well  documented.  Sections  II  and  V 
of  the  current  report  contain  suitable  references  and  examples. 

Not  withstanding  the  development  and  applications  of  metallic 
data  bases  the  growth  of  similar  facilities  for  ceramic  and 
oxide  systems  has  proceeded  at  a  much  slower  pace.  This  was 
certainly  the  case  in  1975  when  the  present  program  was  initiated 
(see  page  6) .  Accordingly,  the  current  work  was  initiated  in 
order  begin  the  development  of  a  data  base  for  oxide  and 
oxynitride  ceramic  systems  which  could  be  applied  to  computation 
of  phase  equilibria  and  stability  over  a  wide  range  of  temperature, 
pressure  and  composition.  One  of  the  objectives  of  this  study 
was  to  develop  the  capability  of  dealing  with  SIALON  systems 
which  are  based  on  combinations  of  Si^N^,  AIN,  A1203  and  Si02. 
Reference  to  Sections  III  and  VI  shows  that  this  objective  has 
been  achieved  through  our  own  efforts  and  those 
which  were  carried  on  by  other  workers  5  l  the  "CALPHAD  group" 

(i.e.  references  2,3  page  42  and  references  3  and  4  page  61). 
Section  4  provides  a  number  of  examples  of  computed  ternary 
sections  through  systems  composed  of  oxynitride  ceramics.  The 
data'Hase "which  has  been  developed  for  the  purpose  of  performing 
the  calculations  displayed  in  Section  4  can  be  employed  to 
•explore  stable  and  metastable  equilibrium  in  multi  component 
oxynitride  systems. 

'  The  three  technical  papers  contained  in  Sections 
II-IV  trace  the  development  of  the  present  data  base  which 
currently  covers  Cr203,  A1203,  MgO,  FeO,  Fe203,  Fe30^,  Si02, 

CaO,  -SijN^  and  AIN.  Section  II  discusses  the  initial  segment 
dealing  with  combinations  of  Cr203,  MgO,  A120j,  FeO  and 

Fe^O^.  The  twelve  quasibinary  systems  and  three  quasiternary 
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systems  were  calculated  and  compared  with  limited  experimental 
data.  The  development  of  lattice  stability,  solution  and  compound 
phase  parameters  for  the  liquid,  spinel,  corundum  and  periclase 
phases  can  be  employed  to  calculate  phase  equilibria  activity, 
vapor  pressure  and  free  energy  of  formation  of  solution  phases 
and  compounds  in  quasibinary  and  quasiternary  systems  composed 
of  the  above  listed  compounds.  The  ternary  isothermal  sections 
which  were  calculated  in  this  work  include  the  Fe^O^-Cr^^-A^O^ 
system  between  1600  and  2400K,  the  MgO-FejO^-A^O^  system  between 
1600K  and  2850°K  and  the  MgO-C^O^-A^Oj  between  1300°K  and 
2800K.  In  the  last  case  the  computed  results  were  found  to 
compare  favorably  with  experimental  results. 

The  second  phase  of  the  development  covered  in  Section  III 
deals  with  inclusion  of  CaO  and  SiC>2.  Each  of  the  quasibinary 
systems  between  their  compounds  and  the  forgoing  group  (i.e.  C^O^, 
A^Oj,  MgO,  FeO,  Fe2Oj  and  Fe^O^)  was  analysed  permitting  calcu¬ 
lation  of  phase  diagram  and  thermochemical  free  energy  of  formation 
of  quasibinary  compounds.  This  description  covered  the  liquid, 
spinel,  periclase,  corundum,  crystobalite,  tridymite  and  the 
alpha  and  beta  quartz  structures  which  occur  in  these  quasibinary 
systems.  In  addition,  ternary  isothermal  sections  were  computed 
between  1600  and  2800K  for  the  MgO-FejO^-Cr^j  system. 

The  extension  to  oxynitride  ceramic  systems  described 
in  Section  IV  was  carried  out  by  developing  lattice  stability, 
solution  and  compound  phase  parameters  for  each  of  the  quasibinary 
systems  between  Si^N^  and  AIN  and  MgO,  A^Oj  and  Si02.  The  beta 
SijN^  and  hexagonal  AIN  structures  were  included  along  with 
the  liquid  and  previously  listed  solid  phase  structures.  The 
melting  point  and  decomposition  temperatures  of  Si^N^  and  AIN 
at  one  atmosphere  pressure  were  calculated  in  the  course  of 
tis  analysis.  Isothermal  sections  in  the  Si02'Al20j-SijN^  and 
Si3N4-Al203-AlN  systems  were  computed  which  are  in  good  agreement 
with  experimental  findings.  Moreover,  the  free  energy  of 
formation  of  the  ternary  compound  phase  X  which  forms  in  the 


SiC^-Al^O^-Si and  the  beta  sialon  phase  of  the  Si^N^-Al^O^-AlN 
system  wore  calculated.  Comparison  of  the  latter  results  with 
independent  findings  obtained  by  the  Stuttgart  group  yields 
good  results  which  can  be  further  improved  by  iteration. 

Section  V  provides  a  listing  of  recent  cased  in  which 
thermochemical  data  have  been  employed  successfully  in  dealing 
with  practical  problems.  Such  documentation  is  difficult  to 
obtain  because  the  secrecy  and  confidentiality  restrictions 
which  surround  economically  viable  processes  necessarily  limit 
publication  of  such  information.  Nevertheless,  such  examples 
come  to  light  over  a  period  of  years.  The  discussion  in  Section 
V  covers  processing  of  in-situ  eutectic  composite  formation 
in  columbium,  nickel  and  cobalt  base  superalloys,  sigma  phase 
formation  in  high  temperature  alloys,  multicomponent  semiconductors 
for  electroptical  devices,  SIALON  ceramics,  molten  salt,  sulfide 
and  oxide  slags  for  ore  reduction  and  refining,  high  temperature 
battery  electrodes,  hardenability  calculations  in  alloy  steels, 
synthesis  of  hydrc^n  storage  alloy  compounds,  grain  refinement 
in  aluminum  alloys,  metatectic  reactions  in  hafnium  alloys  and 
the  growth  of  HEM  silicon  solar  cell  materials. 

The  final  section  of  this  report  discusses  application 
of  the  present  computational  methods  to  the  CVD  processing  of 
doped  germanium-borosilicate  glasses.  Section  VI  illustrates 
the  utilization  of  Data  Bank  information  in  computing  the 
temperature/pressure  range  over  which  the  halides  of  potential 
dopents  can  be  employed  effectively  in  the  CVD  processing  of 
preforms.  As  the  first  step  in  an  analysis  of  these  glasses. 
Section  VI  contains  an  analysis  of  the  Si02-Ge02,  Si02-B203 
and  B203-Ge02  quasibinary  systems  as  well  as  the  Si02-Ge02-B203 
quasiternary  system  over  the  temperature  range  from  1000K  to 
1600K.  Future  activities  will  include  consideration  of  the 
addition  of  Sb20g,  P203  an<*  Ce02  to  system« 


Section  la.  Program  Personnel 

Technical  activities  on  this  program  have  been  performed 
by  L.  Kaufman,  H.  Nesor,  E.P.  Warekois,  D.  Pendse,  J.  Smith, 

G.  Zilberstein,  C.  Thielen,  K.  Taylor,  A.P.  Miodownik,  and 
J.  Nell  during  the  course  of  this  study. 

Section  lb.  Interaction  with  Other  Technical  Activities 

During  the  course  of  this  study  technical  presentations 
covering  this  work  were  given  at  the  ASM/AIME  Materials  Congresses 
in  Chicago,  Illinois  (October  1978)  and  Milwaukee,  Wisconsin 
(September  1979),  CALPHAD  VI-Harwell,  Oxford  England  April  1977, 
CALPHAD  VTI-Stuttgart ,  West  Germany  (April  1978),  Technion- 
Haifa  Israel  (May  1978)  ,  CALPHAD  VUI-Stockholm  Sweden  (May  1979) 
and  National  Physical  Laboratory/Chemical  Society  of  London 
Symposium  on  Industrial  Use  of  Thermochemical  Data,  University 
of  Surrey  Guildford  Surrey  England  (September  1979).  Specific 
interactions  with  scientific  personnel  at  the  Max  Plank  Institute- 
Stuttgart,  West  Germany,  University  of  Surrey,  Guildford  Surrey, 
England  and  at  the  Israel  Institute  of  Technology-Technion  City 
Haifa  Israel  were  most  helpful. 
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Section  II, 


Calculation  of  QuasiBinary  and  QuasiTernary 
Systems  -  I,  CALPHAD  (1978)  2  pp  35-53 


CALPHAD  Vol.  2,  No.  1,  pp.  35-53.  Pergamon  Press,  1978.  Printed  in  Great  Britain. 


CALCULATION  OF  QUASI BINARY  AND  QUASITERNARY 
OXIDE  SYSTEMS  -  I* 


Larry  Kaufman  and  Harvey  Nesor 
ManLabs,  Inc.,  21  Erie  Street 
Cambridge,  Massachusetts  02139,  USA 


ABSTRACT.  A  data  base  is  being  developed  for  calculation  of  quasi-binary  and 
quasi-ternary  phase  diagrams  of  ceramic  systems.  The  initial  seg¬ 
ment  of  this  base  covers  combinations  of  Cr2  03  ,  MgO,  At203  ,  Fe203, 
FejO*and  "FeO".  Lattice  Stability,  Solution  and  Compound  Phase 
Parameters  are  presented  covering  the  liquid,  spinel,  periclase  and 
corundum  phases  in  systems  composed  of  these  oxides.  These  para¬ 
meters  permit  calculation  of  the  thermochemical  properties  and  phase 
diagrams  of  combinations  of  these  oxides  in  good  agreement  with 
limited  experimental  data.  Twelve  calculated  quasi-binary  phase 
diagrams  and  three  calculated  quasi-ternary  phase  diagrams  are 
presented  illustrating  the  applicability  of  the  data  base. 

1.  Introduction 

The  utility  of  computer  based  methods  for  coupling  phase  diagrams  and 
thermochemical  data  for  metallic  systems  has  been  demonstrated  in  many  papers 
published  in  this  journal.  A  considerable  effort  is  being  applied  toward 
development  of  a  suitable  data  base  for  metallic  systems  (1).  Numerous  exam¬ 
ples  have  been  published  showing  how  such  data  on  metallic  systems  can  be 
synthesized  to  successfully  predict  the  thermochemical  properties  and  phase 
relations  in  ternary  metallic  systems  (1-12).  Such  predictions  can  be  applied 
practically  to  explore  conditions  for  in-situ  eutectic  growth,  metallic  glass 
formation,  spinodal  decomposition,  metastable  equilibrium,  and  diffusionless 
transformations.  In  a  recent  workshop  on  "Applications  of  Phase  Diagrams  in 
Metallurgy  and  Ceramics"  held  in  January,  1977  at  the  National  Bureau  of 
Standards  in  Gaithersburg,  Maryland  (13),  many  informational  gaps  in  ceramic 
phase  diagram  data  were  noted  and  the  need  to  develop  computational  techniques 
to  bridge  these  gaps  became  evident.  Such  efforts  have  been  proceeding  indivi¬ 
dually  for  some  time  (14,15).  However,  in  view  of  the  success  achieved  in 
developing  a  tractable  data  base  for  metallic  systems,  it  seemed  to  be  an 
auspicious  time  for  launching  the  development  of  a  similar  data  base  for  ceram¬ 
ic  systems.  This  is  particularly  true,  since  the  need  for  such  information  is 
apparent  (i.e.,  metastable  reactions  are  quite  common  in  ceramic  systems, pres¬ 
sure  is  an  important  variable,  etc.)  and  experimental  data  is  limited  and 
difficult  to  obtain.  Accordingly,  the  present  work  should  be  considered  as  a 
beginning,  much  like  the  initial  attempt  at  estimating  the  lattice  stability 
of  metals  (16,17).  The  only  difference  is  the  experience  gained  with  develop¬ 
ment  of  a  data  base  for  metals  over  the  past  ten  years  and  the  availability  of 
suitable  computational  programs. 


♦This  work  has  been  sponsored  by  the  Air  Force  Office  of  Scientific  Research, 
Bolling  AFB,  Washington,  D.C.,  under  Contract  F44620-76-C-0060. 
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TABLE  1 

Lattice  Stability,  Solution  Phase  and  Compound 
Parameters  Employed  for  Calculation  of  Several 
Metal  -  Oxygen  Systems 

(All  units  in  Joules  per  gra'm  atom  (mole  of  atoms),  T  in  Kelvins) 
Lattice  Stability  Parameters 

FEFEVL  =  371200  -  118. 4T,  ALALVL  =  3012S0  -  107. 9T 

CRCRVL  =  367360  -  120. ST,  OOOOVL  =  3410  -  37. 9T 

MGMGVL  =  133890  -  98. 3T,  OOOOLB  =  230.  -  4.2T 


Solution  Phase  Parameters 


0  to  50  a/o  Oxygen 


SO  to  100  /o  Oxygen 


LFEOO 

0 

423420 

- 

4S6.1T 

LOOFE 

-1000810 

+ 

12S.ST 

-1424230 

- 

S81.6T 

LCROO 

0 

-882820 

+ 

117. 2T 

L00CR 

-176S6S0 

■f 

234. 3T 

-882820 

+ 

117. 2T 

LMGOO 

-90790 

• 

411 . 3T 

-90790 

- 

411. 3T 

LOOMG 

-1940120 

+ 

626. 3T 

-1940120 

+ 

626. 3T 

LALOO 

-1287830 

+ 

147. 7T 

-1287830 

+ 

147. 7T 

LOOAL 

-1287830 

+ 

147. 7T 

-1287830 

+ 

147. 7T 

for  each  solution  phase  with  u;  ic|um.tu  \.y  uj  ,  uj  m.., 

E'hcp  phases  etc.,  and  with  LIJ^  replaced  by  BJJ  and  EIJ  etc.  In  the  present 
symbolism,  one  uses  the  symbol  FEFELB  to  designate  the  free  energy  of  pure 
liquid  iron  minus  the  free  energy  of  pure  bcc  form  of  iron.  Thus  FEFELB=  -FEFEBL 
Similarly,  the  free  energy  of  pure  liquid  chromium  minus  the  free  energy  of 
pure  bcc  chromium  is  designated  by  CRCRLB.  If  I=iron  and  J=nickel  in  Equation 
(1),  then  Table  10  of  Reference  1  shows  that  LFENI=  -8368  +  2.72T  (J/g.at.)  and 
LNIFE=  -32217  +  9.205T  (J/g.at.).  Figure  1  and  Appendix  I  shows  how  the  free 
energy  of  a  solution  phase  can  be  described  along  a  specific  quasi-binary  join 
in  a  given  ternary  system.  This  approximation  is  used  here  to  reduce  the  ter¬ 
nary  description  to  a  quasi-binary  description  in  order  to  employ  the  informa¬ 
tion  stored  in  the  data  bank  to  deal  with  high  order  systems.  The  free  energy, 
G’f’,  of  a  compound  phase  i.,  J  is  defined  by  Equation  (2),  (see  Equation  .14 
page  60  of  (1))  as  x* 

6*  =  (1-xJ  °G®  +  xt  *G  ?  +  x.Cl-x,)  [(1-xJLIJ  i  x.LJI-C'Tj] 

where  6  refers  to  the  base  phase  selected  for  the  compound  and  C IT]  (which  is  a 
temperature  dependent  function)  is  the  compound  parameters.  Taoles  1;4  define 
the  iron-oxygen,  chromium-oxygen,  aluminum-oxygen,  magnesium-oxygen,  iron- 
magnesium,  chromium-magnesium,  iron- aluminum,  iron-chromium,  alummum-chromiur. 


•G-*  replaced  by  °Gi 


yl 


etc. ,  for  the  B=bcc, 
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Compound 

(Base  Phase  is  BCC) 

Structure 

Compound 

Parameter 

Fe0.S°0.S 

Periclase  (P) 

28030  + 

12.S5T 

1 

Fe0 . 429°0 . 571 

Spinel  (S) 

40170  + 

11.63T 

| 

Fe0.4°0.6 

Corundum  (C) 

51040  + 

11.71T 

1 

Cr0.4°0.6 

Corundum  (C) 

112130  - 

8.37T 

f 

M^0 . S°0 . 5 

Periclase  (P) 

182004  - 

16.74T 

[ 

^^0 .4^0.6 

Corundum  (C) 

96650  - 

2.51T 

\ 

1 

i 
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TABLE  2 


Comparison  of  Calculated 

and  Observed  (18,20)  Free 

Compound 

Energy  of  Formation  of  Binary  Oxide  Compounds 

(All  units  in  Joules  per  mole,  T  in  Kelvins) 

(P*Periclase ,  S»Spinel,  C'Corundum) 

Free  Energy  of  Formation 

FeO  (P) 

Calculated 

-267860  +  67. 1ST 

Observed 

-271120  +  66.94T 

Fe304  (S) 

-1166600  +  382. 21T 

-1086S80  +  292. 88T 

Fe203  (C) 

-894330  +  322. 59T 

-803328  +  251. 04T 

ft 

O.SFeO+O.lZSFe  +  0.12S  Fe,0.;  AG  - 

-11900  +  14.18T;  TQ  *  839K 

0.133Fe304^-0. 

03340^^0 . 2Fe20j ;  AG  B 

-23700  +  13.71T;  TQ  *  1726K 

Cr203  (C) 

Calculated 

-1204990  +  276. 77T 

Observed 

-1129680  +  251. 04T 

T*1000 

-874460 

-928220 

1500 

-747260 

-789830 

2000 

-618810 

-651450 

Calculated 

Observed 

ai2o3  (C) 

T*i000 

-1359800 

-1359800 

1500 

-1202900 

-1196620 

2000 

-1046000 

-1035540 

Calculated 

Observed 

MgO  (P) 

T»  500 

-552960 

-547520 

1000 

-499360 

-492460 

1500 

-426810 

-423840 

2000 

-322210 

-319240 

2500 

-217610 

-215060 

3000 

-113010 

-110460 

TABLE  3 

Solution  Phase  Parameters  for  Several  Binary 

Metal  Systems 

(All  units  in  Joules  per  gram  atom  (mole  of  atoms),  T  in  Kelvins) 

LFEHG 

C 

64850 

LFECR 

87000  -  10.46T 

LMGFE 

M 

64850 

LCRFE 

= 

87000  -  10.46T 

LCRMG 

St 

5'1050 

LALCR 

£ 

-46440 

LMGCR 

e 

51050 

LCRAL 

as 

-46440 

LFEAL 

* 

-62760  -  4.60T 

LALMG 

= 

-10040  +  3.35T 

LALFE 

= 

-96230  +  33.05T 

LMGAL 

St 

-10040  +  3.3ST 
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TABLE  4 

Ternary  Liquid  Solution  Phase  Parameters 
(For  Several  Metal^Metal-Oxygen  Systems) 

(All  units  in  Joules  per  gram  atom  (mole  of  atoms,  T  in  Kelvins) 


System 

TRNL . 

System 

TRNL 

AE-O-Cr 

0 

Fe-O-AE 

1757280  -  1125. ST 

Fe-O-Cr 

1439300  -  912. IT 

AE-O-Mg 

1882800  -  1121. 3T 

Cr-O-Mg 

2217S20  -  12SS.0T 

TABLE  5 

Lattice  Stability 

Parameters 

for  Oxides 

(All  units  in 

Joules  per  gram 

atom  (mole  of  atoms),  T  in  Kelvins) 

Symbol  Name 

M 

Example 

C 

Corundum 

HR  10 

A*2°3 

S 

Spinel 

CF  56 

Fe3°4  CCo3°4) 

P 

Periclase 

CF  10 

MgO  (NaCE) 

Abbreviations 

W0  •  Fe0.S°0. 

5>  FO  *  FeQ 

.4°0.6'  H0  *  Fe0.429°0.571 

AO  «  AEq^Oq^ 

6*  C0  *  Cr0 

. 4°0 . 6 ’  M0  ’  M^0 . 5°0 . 5 

WOWOLC  - 

7280 

-  7.28T 

AOAOLC 

-  23640  -  10.21T 

NOWOLS  - 

9017 

-  6.44T 

AuAOLS 

«  19205  -  9.08T 

KOWOLP  « 

14016 

-  8.54T 

AOAOLP 

-  4477  -  12.84T 

FOFOLC  ■ 

18117 

-  8.37T 

COCOLC 

*  343S0  -  13.47T 

FOFOLS  » 

153SS 

-  7.29T 

COCOLS 

»  22594  -  12.S5T 

FOFOLP  - 

9414 

-  6.28T 

COCOLP 

*  34058  -  15.48T 

HOHOLC  - 

9707 

-  9.71T 

MOMOLC 

-  22594  -  12.55T 

HOHOLS  - 

15816 

-  8.45T 

MOMOLS  . 

.«  19560  -  10.46T 

H0H0LP  - 

12217 

-  7.20T 

MOMOLP 

-  47698  -  15.40T 

and  aluminum-magnesium  systems  in  order  to  apply  the  description  shown  in 
Figure  1  and  Appendix  I.  In  these  Tables,  L  refers  to  the  liquid  phase  while 
B  and  V  refer  to  the  bcc  and  vapor  phase  respectively.  The  symbol  MG  is  used 
for  magnesium,  AL  for  aluminum  and  00  for  l/202.  Thus,  the  basis  in  all  cases 
is  one  mole  of  atoms  or  one  gram  atom.  The  numerical  values  chosen  for  the 
lattice  stability,  solution  and  compound  phase  parameters  were  derived  from 
relevant  thermochemical  and  phase  diagram  data  for  oxygen-containing  systems 
(18-24)  supplemented  by  previous  assessments  of  lattice  stability  (1,17)  and 
solution  phase  parameters  (i.e.,  iron-chromium  (1,2),  iron-aluminum  and 
aluminum-chromium  (2),  aluminum-magnesium  (9)).  The  selection  of  the  ternary 
liquid  interaction  parameters  for  the  aluminum-oxygen-chromium,  iron-oxygen- 
chromium,  chromium-oxygen-magnesium,  iron-oxygen-aluminum  and  aluminum-oxygen- 
magnesium  systems  shown  in  Table  4  was  made  on  the  basis  of  the  data  contained 
in  Reference  (20).  The  description  of  Tables  1-4  permits  calculation  of  the 
thermochcmical  properties  of  the  above  noted  oxide  systems  and  calculation  of 
their  phase  diagrams  under  standard  (one  atom  phase  pressure)  conditions.  In 
addition,  it  permits  calculation  of  the  excess  free  energy  of  the  liquid  phase 
in  the  quasi-binary  join  of  a  ternary  metal-oxygen-metal  system  as  shown  in 
Figure  1  and  Appendix  I.  This  facility  will  be  employed  subsequently.  The 
calculated  iron-oxygen  and  chromium-oxygen  phase  diagrams  resulting  from  this 
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Figure  5.  Calculated  Al0.4<>0.6'Fe0.4O0.6  Phase  Diagram. 

(Note:  Figure  2  shows  thet  FeQ  ^Oq  6  Decomposes 
into  Fe^  429°o  571  and  #I)0VC  l726*0- 


description  are  shown  in  Figures  2  and  3.  These  diagrams  are  in  reasonably  . 
good  agreement  with  the  observed  diagrams  if  one  allows  foT  the  idealizations 
inherent  in  treating  iron  oxide  (wustite)  as  a  line  compound  at  50/50  stoichio¬ 
metry. 

3.  The  Lattice  Stability,  Solution  and  Compound  Phase  Parameters 
for  the  Spinel,  Periclase  and  Corundum  Structures 

Table  5  shows  the  values  derived  for  the  lattice  stability  of  the  liquid, 
L,  spinel,  S,  periclase,  P,  and  corundum,  C,  structures  of  the  oxides  AJt303, 
Cr20j,  MgO,  Fe30,,,  Fe203  and  "FeO".  These  lattice  stability  values  were 
employed  to  compute  the  quasi-binary  phase  diagrams  shown  in  Figures  4-15.  It 
should  be  noted  that  the  lattice  stability  values  are  given  for  a  gram  atom 
(one  mole  of  atoms)  of  each  of  the  oxides!  Moreover,  examination  of  Table  5 
shows  that  a  gram  atom  of  each  of  the  oxides  (i.e.,  Cr0.uO(,.6)  has  been  assign¬ 
ed  a  specific  abbreviation  (i.e.)  CO  *  Cro.nOo.6.  The  g-atom  free  energy  of  the 
liquid  minus  the  corundum  form  of  Cro.nOj.g  is  designated  by  the  symbol  COCOLC 
and  equals  34350  -  13.47T  Joules  as  shown  in  Table  5.  In  direct  analogy  to  the 
discussion  of  solution  phase  parameters  in  Section  2  above,  the  solution  phase 
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TABLE  6 

Quasibinary  Solution  Parameters  for  Oxide  Systems 
(All  units  in  Joules  per  gram  atom  (mole  of  atoms) ,  T  in  Kelvins) 


LAOCO  = 

-7422  + 

13.05T 

LAOFO  = 

-14297  + 

5.6ST 

SAOCO  * 

-7422  + 

13.05T 

SAOFO  = 

-14297  + 

5.6ST 

PAOCO  = 

-7422  + 

13.05T 

PAOFO  = 

6623  + 

S.6ST 

CAOCO  = 

-7422  + 

13/05T 

CAOFO  = 

6623  + 

S.65T 

LAOHO  = 

-203S0 

+ 

11.2ST 

LAOMO  = 

-33530  + 

22.93T 

SAOHO  = 

-20350 

+ 

11.25T 

SAOMO  = 

-33S30  + 

22.93T 

PAOHO  = 

570 

+ 

11.25T 

PAOMO  = 

29230  - 

4.69T 

CAOHO  = 

570 

+ 

11.25T 

CAOMO  = 

-12610,  + 

22.93T 

LAOWO  = 

-39330 

+ 

2S.73T 

LCOAO  = 

-7422  + 

13.05T 

SAOWO  = 

-39330 

+ 

2S.73T 

SCOAO  = 

-7422  + 

13.0ST 

PAOWO  = 

-18410 

+ 

25.73T 

PCOAO  = 

-7422  + 

13.0ST 

CAOWO  = 

-18410 

+ 

25.73T 

CCOAO  = 

-7422  + 

13.05T 

LCOFO  - 

-26080 

+ 

19.16T 

LCOHO  = 

-37780  + 

26.19T 

SCOFO  = 

-26080 

+ 

19.16T 

SCOHO  = 

-37780  + 

26.19T 

PCOFO  = 

-26080 

+ 

19.16T 

PCOHO  = 

-4310  + 

15.73T 

CCOFO  = 

7390 

+ 

8.70T 

CCOHO  = 

-4310  + 

15.73T 

LCOMO  = 

780 

+ 

7.45 

LCOWO  = 

-70760  + 

44.22T 

SCOMO  = 

780 

+ 

7.45 

SCOWO  = 

-70760  + 

44.22T 

PCOMO  = 

21700 

+ 

3.26 

PCOWO  = 

-49840  + 

44.22T 

CCOMO  = 

21700 

+ 

7.45 

CCOMO  = 

-37290  + 

33.76T 

LFOAO  = 

58150 

. 

30.03T 

LFOCO  = 

41000  - 

18.49T 

SFOAO  = 

58150 

- 

30.03T 

SFOCO  *= 

41000  - 

18.49T 

PFOAO  = 

79070 

- 

30.03T 

PFOCO  = 

41000  - 

18.49T 

CFOAO  = 

79070 

30.03T 

CFOCO  = 

41000  - 

18.49T 

LFOMO  = 

-513S0 

+ 

26.82T 

LHOAO  = 

59390  - 

35.77T 

SFOMO  « 

-£1350 

+ 

26.82T 

SHOAO  = 

59390  - 

35.77T 

PFOMO  = 

-9510 

+ 

14.27T 

PHOAO  * 

80310  - 

35.77T 

CFOMO  = 

-30430 

+ 

26.82T 

CHOAO  = 

80310  - 

35.77T 

LHOCO  = 

37810 

. 

15.44T 

LHOMO  = 

-30040  + 

19.75T 

SHOCO  = 

37810 

- 

1S.44T 

SHOMO  = 

-34220  + 

19.75T 

PHOCO  = 

37810 

- 

15.44T 

PHOMO  = 

6360  + 

11.38T 

CHOCO  = 

37810 

* 

15.44T 

CHOMO  = 

-9120  + 

19.75T 

LMOAO  « 

59910 

- 

27.11T 

LMOCO  = 

87520  - 

39.96T 

SMOAO  = 

59910 

- 

27.11T 

SMOCO  = 

87520  - 

39.96T 

PMOAO  = 

122670 

- 

54.73T 

PMOCO  = 

108440  - 

44.14T 

CMOAO  » 

80830 

- 

27.11T 

CMOCO  = 

87520  - 

39.96T 

LMOFO  » 

-22340 

+ 

13.56T 

LMOHO  = 

-10370  + 

10.54T 

SMOFO  = 

-22340 

+ 

13.56T 

SMOHO  = 

-14556  + 

10.54T 

PMOFO  « 

18500 

+ 

1.00T 

PMOHO  = 

26030  + 

2.18T 

CMOFO  = 

-2420 

+ 

13.56T 

CMOHO  =• 

10550  + 

10.54T 

LMOWO  - 

16890 

- 

4.06T 

LWOAO  == 

63180  - 

31.71T 

SMOKO  = 

16890 

- 

4.06T 

SWOAO  = 

63180  - 

31.71T 

PMOWO  = 

8520 

- 

4.06T 

PWOAO  = 

84100  - 

31.71T 

CMOWO  «= 

8520 

- 

4.06T 

CWOAO  = 

84100  - 

31.71T 
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TABLE  6  (continued) 


LWOCO  -  31440  -  10.00T 
SWOCO  -  31440  -  10.00T 
PWOCO  -  50840  -  10.00T 
CWOCO  -  31440  -  10.00T 


LWOMO  -  16800  +  3.93T 
SWOMO  -  16800  +  3.93T 
PWOMO  -  8430  +  3.93T 
CWOMO  -  8430  +  3. 93T 


TABLE  7 

Quasibinary  Compound  Parameters  for  Oxide  Systems 
(All  units  in  Joules  per  gram  atom  (mole  of  atoms)  T  in  Kelvins) 


Compound 

Structure 

Base 

Compound  Parameters 

(M0.286°0.428) 

*(Mg0.143°0.143) 

(SP) 

(S) 

74050 

(AA0.286°0.428) 

4(Fe0.146°0.143^ 

(SP) 

(S) 

40580  + 

7.11T 

(Cr0.286°0.428) 

•(Mg0.143°0.143> 

(SP) 

(S) 

134720  - 

7.11T 

(Cr0.286°0.428:) 

'(Mg0.143°0.143) 

(SP) 

(S) 

38910  ♦ 

16.74T 

(Fe0.286°0.428) 

'(Mg0.143°0.1435 

(SP) 

(S) 

18410  + 

13.72T 

^VsoeS^os5 

•(Mg0.143°0.143) 

(SP) 

(S) 

O 

T*BLE  8 

Comparison  of  Calculated  and  Observed  (2(5  Free  Energies 
of  Formation  of  Complex  Spinel  Oxides 

(P*Periclase,  C*Corundum,  SP»Spinel) 

Reaction  Free  Energy  of  Formation.  Joules ,  T  in  Kelvins 


FeO(P) 

+  Ai,0T(C)-Fe0.A£,<X.(SP) 

-40070  - 

6.74T 

(CALC) 

4  0  4  0 

-39740  - 

7.12T 

(OBS) 

MgO(P) 

+  Cr,0,(C)-*-Mg0'Cr,0.  (SP) 

-40770  - 

13.03T 

(CALC) 

4  0  4  0 

-41120  - 

14.00T 

(OBS) 

MgO(P) 

+  A*,0t(C)-MgO-At,O,(SP) 

-36961  - 

3.22T 

(CALC) 

4  0  4  0 

-35560  - 

2.05T 

(OBS) 

FeO(P) 

+  Cr,0_(C)-*’Fe0*Cr,0,(SP) 

-46100  + 

8.32T 

(CALC) 

4  0  4  0 

-45190  + 

8.38T 

(OBS) 

MgO(P) 

+  Fe,0.(C)+Mg0*Fe,0-(C) 

-20940  - 

2.20T 

(CALC) 

4  0  4  0 

-19330  - 

2.0ST 

(OBS) 

parameters  for  one  gram  atom  of  the  liquid  on  the  AtjOj-MgO  join  of  the 
AA-O-Mg  system  are  designated  by  LAOMO  and  LMOAO .  These  parameters  have  been 
computed  from  Tables  1-4  and  the  procedure  outlined  in  Appendix  I,  and  are  also 
displayed  in  Table  6.  Reference  to  Table  6  shows  that  LAOMO  •  -33S30  +  22.93T 
and  LMOAO  -  59910  -  27.11T  respectively.  ' 

Tables  6-8  summarize  the  solution  phase  parameters  for  the  remaining 
phases  (i.e.,  C,  corundum,  S,  spinel,  P,  periclase).  Thus  the  solution  para¬ 
meters  for  one  gram  atom  of  the  corundum,  solution  in  the  AAj.»0.t-Cro.<,00.«  j°in 
of  the  Ai.-0-Cr  system,  can  be  read  as  CAOCO  ■  *7422  ♦  13.Q5T  and  CCOAO  *  -7422  + 
13. 0ST  Joules  respectively.  Tables  7  and  8  display  the  compound  parameters 
employed  in  calculating  the  quasi-binary  phase  diagrams  shown  in  Figures  4-15. 
Table  8  shows  that  these  same  parameters  can  be  employed  to  generate  the  free 
energy  of  formation  of  complex  spinel  oxides  in  agreement  with  observation  (20). 
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Thus  the  definitions  of  Equations  (1)  and  (2)  coupled  with  the  quasi-binary 
approximation  shown  in  Appendix  I  and  parameters  shown  in  Tables  1-8  generate, 
the  one  atmosphere,  thermochemical  properties  and  phase  diagrams  for  the 
quasi-binary  joins  displayed  in  Figures  4 - 1 S . 

4 .  Discussion  Of  The  Results  Of  The 
Quasi-Binary  Approximation: 

Table  5  provides  the  "liquid-solid"  free  energy  differences  (and  the 
lattice  stability  parameters)  for  the  six  oxides  of  current  interest.  Thus, 
AOAOLC  is  equal  to  the  free  energy  of  liquid  At,  ..O,  s (i . e. ,1 /5A120,)  minus 
the  free  energy  of  the  corundum  form  of  AJt0 .  „0, | 1 .  *  Information  on  the  heat 
of  fusion  for  the  stable  forms  of  each  oxide  was  taken  from  review  publications 
and  recent  experimental  studies  (18,20,25).  As  indicated  earlier,  AOAOLS  and 
AOAOLP  represent  the  corresponding  free  energy  differences  between  the  liquid 
and  spinel  and  the  liquid  and  periclase  forms.  This  description  permits  cal¬ 
culation  of  the  free  energy  differences  between  the  C,  S  and  P  forms  since 
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the  free  energy  of  the  spinel  form  minus  the  free  energy  of  the  corundum  form 
of  Af.0  ..Oo  6>  AOAOSC,  is  equal  to  AOAOLC  -  AOAOLS  =  4435  -  1.13T  J/g.at. 

The  enthalpy  difference  of  4435  J/g.at.  is  in  keeping  with  experimentally 
observed  results  (26,28).  A  similar  conclusion  can  be  drawn  for  the  enthalpy 
difference  of  2762  J/g.at.  between  the  spinel  and  corundum  forms  of  Feo.tOo.s 
derived  from  the  difference  FOFOLC  -  FOFOLS  given  in  Table  5  (26,28).  The 
lattice  stability  parameters  shown  in  Table  5  were  derived  by  using  experimen¬ 
tal  values  as  a  base  and  "filling  in"  t’ s  missing  components  by  performing 
an  "ideal  solution"  iteration  on  Figure:.  4-15  as  in  the  case  of  metals  (16,17). 
The  solution  phase  parameters  for  the  liquid  shown  in  Table  6  are  fixed  by 
the  procedure  outlined  in  Appendix  I  and  the  results  for  the  metal- oxygen, 
metal-metal  and  metal-oxygen-metal  systems  shown  in  Tables  1-4.  Reference  to 
Table  6  discloses  that  the  corresponding  values  for  the  S,  P  and  C  solid 
solutions  are  usually  equal  to,  or  more  positive  than  the  L  parameters.  The 
S,  spinel  solution  parameters  are  usually  equal  to  the  liquid  parameters  while 
the  periclase,  P,  and  corundum,  C,  parameters  are  usually  more  positive  than 
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L  and  S  (and  equal).  These  smal] 
differences  were  obtained  by 
iterative  calculation  of  Figures 
4-15  starting  with  equal  values 
of  L,  S,  P  and  C. 

Finally,  the  quasi-binary 
compound  parameters  shown  in 
Table  7  were  fixed  to  yield  con¬ 
formity  with  the  phase  diagram 
and  the  observed  thermochemical 
data  shown  in  Table  8. 

5.  Calculation  of  the 
Activity  of  Quasibinary 
and  Ternary 

Components  of  Solution  Phases 

.  Defining  the  activity  of 
quasi-binary  components  across  a 
given  join  is  quite  straight¬ 
forward  (2)  based  on  the  defini¬ 
tions  of  Equation  (1).  Thus, 
the  activity  of  CO  and  AO  in  the 
C  phase  of  the  CO-AO  system 
(Fig.  4)  where  x  is  the  fraction 
of  AO  and  the  reference  states 
are  pure  corundum  CO  and  pure 
corundum  AO  is  given  by  Equa¬ 
tion  (3)  and  (4)  as 

c 

RT*naco  s  RTf.n(l"X)  (3^ 

4x2 [CC0A0*(2x-l) (CA0C0-CC0A0) ] 


RT*naAQ  =  RTfcnx  ^ 

4  (1-x)2  [CCOAO+2x(CAOCO-CCOAO) ] 


/o  MgO  -V 
38  48 


1871 

'  Spinel 
Solid  / 
1000  Solutlo/ 

<S>/ 


O.4:jw0.S71 
Figure  14. 


Periclase 

(P) 


Calculated  Fc0  429<>0  S71-MS()  $O0  5 
Phase  Diagram. 

¥/o  MgO  — 

20  2R  36  46  S7  70  8 


L 

- 

■  1642 

p 

- 

TABLE! 

671 

....  7N 

Figure  15.  Calculated  Feo.S°O.S'M8C.5°O.S  phase 
Diagram.  (Note:  Figure  2  shows  that 
Fe0  j  Decomposes  into  Fe  and 

Fe0.«29°0.S71  *clow 


Similar  expressions  can  be  F  0  Sclow  ,J9K). 

generated  for  each  of  the  phases. 

Derivation  of  equations  for  the 
activity  of  the  ternary  compo¬ 
nents  (i.e.,  Cr,  O  and  A*)  can  only  be  carried  out  for  the  liquid  phase  where 
the  free  energy  has  been  explicitly  stated  in  terms  of  ternary  components  (and 
not  in  terms  of  end  members  on  a  given  binary  join).  To  illustrate  this  proce¬ 
dure,  consider  the  liquid  At-O-Mg  case  where  it  is  desired  to  compute  RTtnaAi, 
RTlna!'.  and  RTtna!(  .  In  the  case  of  the  first  two  components,  a*.  =p*s/0Pa5  and 

I  L  U  L  L  Ax  ax  ax 

aM  =  PM  ^°PM  w^ere  P^l  anc*  pm  are  vaPor  pressures  of  A1  and  Mg  over  the 


liquid  oxide  and  °pAi  and  °pjjg  are  the  vapor  pressures  of  the  pure  metals  at 
the  temperature  in  question.  A  similar  result  holds  for  the  oxygen  with  a 
slight  complication  due  to  the  use  of  ( 1  / 2 ) 02  ass  component!  Thus  a^fpg  ^°p0 
and  since  0.5RTfcn°p  ^  =  OOOOi.V  =  -3410  ♦  37. 9T  at  one  atmosphere,  then 
RTtnp^  =  -6820  +  75. 8T  +  2RTtnag  ,  where  p£  is  the  pressure  in  atmospheres.  The 

equations  for  component  activities  and  vapor  pressures  can  now  be  written  ex¬ 
plicitly  as  Equations  (5) -(7)  for  the  At-O-Mg  case  where  x  is  the  fraction  of 
0  and  y  is  the  fraction  of  Mg.  Note  that  along  the  A0.  4O0. 6'M8o • 5°  0.5  j°in  t^e 
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Bquations°(S) ••  (7)  1appfy°toXalI^*>valuesn0f  x  2nfyX J7) ;6*°*2y*  Nev”theless 


KJ'Zn^-x-y)  *  x  LALOO  +,y2LALMG  +  xy[LALOO  +  LALMG  -  LOOMG; 


"  X  '1]  (UmL-UW°)  *  -1]  (LMGAL- LALMG) 

(1-x) 


^7  CLMGOO-LQ0^  •  ty [l -2 (x+y) ]TRNL(A£-0-Mg) 


*  (1-X)2UL0°  *  y^LALMG-y(l-x)  [LALOO  *  LALMG  -'loOMG] 

♦  X2r*(2-y)  (2-3x-2vO  ,  , 

a-y)2  x(l-y)  "  +  y  fjV  (LMGAL-LALMG) 


<fep[y(l-x)-x  1  (UK00-L00^  +  fy(l-y)-xyC3-2Cx.y))]  TRNL  (AL-O-Mg) 


^  "  m"y  +  x2wL0°  +  2LALMG-x(l-y)  f.ALOO  ♦  LALMG  -  LOOMG] 

+X  fT V"  "1J  (“QM-UWO)  +  y2(Zl2-*)  +  (2-3y-2x) 

/i  --\2  ..  .  *iJ 

.  2 


(LMGAL- LALMG) 


(I-X)  y(l-x) 

2-[2x+y-xy.y2]  ( IMGOO- LOOMG)  +  fx(l-xl  vVrx  ?r  », 

(x>y)  lXU  *xy (3-2  (x+y)  )  ]  TRNL  (AJt-O-Mg) 


(S 


(6 


(7] 


These  equations  car 
“?  emPf?y?d  to  calculate 
the  activity  and  vapor 
pressure  of  each  elemen¬ 
tal  component  in  the 
liquid  phase  at  any  com- 
position.  Conversely 
they  can  be  used  to  cal¬ 
culate  the  composition  a 
a  given  oxygen  pressure.’ 


Calculation  of 
Quasi-Ternary  Phase 
Diagrams 


I 


The  extension  to 
2-a!i*j*I'nary  systems  c 

accomplished 
oy  letting  each  of  the 

W*ntS,Le‘'  C0>  AO 

®?c* »  play  the  role 
of  the  elements  in 
Appendix  I  and  Referenc< 
C2)*  figures 
16-21  show  the  resultini 
sections  com¬ 
puted  for  the  Fe.o*- 

At*n*A*3°£’  M8°'Fe jOv- 

AtiO,  and  Mg0-Cr203- 
AtjO,  at  selected  tem¬ 
peratures.  in  line 
with  previous  practice. 
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the  free  energy  of  ternary  solution  phases  are  synthesized  from  their  binary 
edge  components  by  using  Kohler's  Equation  (Appendix  I)  with  a  zero  ternary 
interaction  parameter.  The  free  energy  of  ternary  compound  phases  is  defined 
along  specific  joins  (see  Reference  2,  p  2117,  Eqs.  (1,8).  As  an  example,  we 
consider  the  spinel  compound  phase  in  the  system  AO-MO-CO  (see  Figure  21). 

This  compound  runs  between  the  At2Q3‘MgO  and  Cr203*Mg0  compositions  designated 
(Table  7)  as  (Ato  •  2  8«Oo .  •» 2 »)  •  (Mgo .  1  <.  3  Oo .  1 4  3)  and  (Cro  •  2  asOo  •  42 »)  ’  (Mgo  •  j  4  aOo  •  1  hi). 
If  x  is  the  fraction  MO  and  y  is  the  fraction  CO  then  the  free  energy  of  the 
spinel  (SP)  compound  AO,.  l.MO  1  -  CO,.  „  .M0„  is  given  by  Equation  (8)  as 


GSP  .  (l-x-r)GP0  .  .  y{^0  •  m[ylny»(l-»-y>ln(l-x-,),.(l-x)uiO-.)] 

*  I'  '  &!“«  *  <*>“»  *  “e 

where 

AG£  «  0,  AGa  *  (l-xi)xl  t(l-xi)  LAOMO  ♦  xlLMOAO  -  CA] 

AGg  *=  (l-xjx,  ((1-xJ  LCOMO  ♦  x,  LMOCO  -  Cg] 

and  ^  1  1 

*  =  x.  +  yp  ;  p  «  (x,  -  x*)/(i-x,),  x;  =  0.286,  x,  *  0.286,  p  =  0 

CA  =  74050  and  Cfi  =  38910  +  16.74T  (From  Table  8) 


(8) 

I 

(9) 

(10) 


7.  Discussion  of 
Results 


The  present  description 
of  the  iron-oxygen  system 
summarized  in  Tables’  1  and  2 
which  generate  Figure  1,  lead 
to  the  calculated  expressions 
for  the  free  energy  of  forma¬ 
tion  of  "FeO",  FeaOi,  and 
Fe20j  shown  in  Table  2.  These 
calculated  free  energies  of 
formation  can  be  employed  to 
examine  the  decomposition  of 
"FeO"  into  iron  and  Fe30*  on 
cooling  or  the  decomposition 
of  FeaOj  into  Fe30i,  and  oxy¬ 
gen  gas  on  heating.  Table  2 
shows  these  calculations 
which  are  the  basis  upon 
which  the  reaction  isotherms 
shown  in  Figure  2  were  con¬ 
structed.  Any  self  consis¬ 
tent  system  of  thermochemical 
data  for  this  system  would  be 
required  to  generate  the 
appropriate  values  of  T0  for 
these  decomposition  reactions 
(i.e.,  839K  and  1726K  respec¬ 
tively).  Nevertheless  in 
many  instances  ,■  this  test  is 
not  made! 

Figure  2  shows  that 
"FeO"  is  not  stable  below 
839K.  Consequently  the 
quasi-binaiy  joins  which 
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include  the  Fe0  5O0>5 
(i.e.,  WO)  component 
(e.g.,  Figs  8,  12  15) 
are  metastable  with 
respect  to  the  decom¬ 
position  of  WO  below 
839K.  In  a  similar 
fashion.  Figure  2 
shows  that  Fe0.»O0.6 
decomposes  above 
1726K.  Thus  each  of 
the  joins  containing 
the  component  (FO)  is 
metastable  above  1726K. 

As  a  consequence, 

Figures  5,  9  and  13 
illustrate  this  meta¬ 
stability.  Thus  the 
lattice  stability, 
solution  and  compound 
phase  parameters  for 
the  elemental  compo¬ 
nents,  and  their 
binary  combinations, 
shown  in  Tables  1-3, 
were  employed  to 
characterize  the 
thermochemical  pro¬ 
perties  and  the  phase 
diagrams  of  the  iron-  Figure  14.  Calculated  Isothemal  Sections  in  the  Mlo.5Vs‘FeO.«90O.S7r 

oxygen,  aluminum-oxygen,  A‘o.«°o.6  Systen‘ 

magnesium -oxygen,  chro¬ 
mium-oxygen,  iron- 

magnesium,  chromium-  magnesium,  iron -aluminum,  i ron- chromium ,  aluminum-chromium 
and  aluminum-magnesium  systems.  These  binary  systems  form  the  edges  of  the 
At-O-Cr,  Fe-O-Cr,  Cr-O-Mg,  Fe-0-A2.  and  AX.-0-Mg  ternary  systems  in  which  the 
quasi-binary  "joins”  of  interest  lie.  Therefore,  the  development  of  the 
present  base  is  built  on  defining  the  three  edge  binary  systems  (i.e.,  genera¬ 
ting  a  small  number  of  lattice  stability,  solution  and  compound  parameters 
capable  of  being  used  to  compute  the  binary  phase  diagrams  and  thermochemical 
properties)  and  then  using  the  quasi-binary  approximation  shown  in  Appendix  I 
to  develop  a  description  of  the  liquid  free  energy  on  a  specific  quasi-binary 
"join"  in  the  ternary!  The  results  of  this  procedure  is  shown  in  Table  6. 

These  results  are  then  combined  with  the  lattice  stability  values  for  the 
quasi-binary  components  shown  in  Table  5  and  the  remaining  solution  phase  para¬ 
meters  shown  in  Table  6  to  compute  the  "quasi-binary"  phase  diagrams  shown  in 
Figures  4-15.  Stable  compound  phases  (mainly  spinels  (SP)  are  added  by 
defining  a  base-phase  (S)  and  a  compound  parameter  as  shown  in  Table  7.  The 
resulting  free  energy  of  the  compound  must  be  consistent  with  the  observed  free 
energy  of  formation  of  the  complex  compound  from  it’s  "quasi-binary"  components 
and  the  quasi-binary  phase  diagram.  The  means  employed  to  define  the  free 
energy  of  the  quasi-binary  compound  is  identical  to  that  employed  for  binary 
compounds  (1-5,  17)  and  can  be  seen  explicitly  by  considering  the  case  of  the 
spinel  AOo . 7i  1.MO0 . 2»6 »  which  can  be  described  as  (1/7)  (Ai^Os-MgO)  or 
(At*.2««0# . *2*) * (Mgo . maOo . 1 m) .  Reference  to  Tables  6  and  7;  Figure  7  and 
Equations  8  and  9  illustrate  how  the  present  definitions  fix  the  free  enerev 
of  this  compound.  For  example,  the  discussion  of  Equations  8-10  deals  with  the 
AOo.7j»MOo.2»«  -  COo.7mMOo.2»s  join  in  which  the  fraction  of  CO  is  equal  to  y. 
When  y  *  0 ,  Equation  (8)  provides  a  description  of  the  free  energy  of  the 
AOo  ,7i 4MO0.  2«s  (SP)  spinel  which  is  equal  to: 

GSP  »  0.714G^q  +  0.286Gj^  +  (0.714) (0.286) [0.714LA0M0  +  0.286LM0AO  -  C] 


(ID 


SP*  *AlO.286°0.42JI'  fC0. 50t»°O.408) 


S«SpineI 
Sol  ill 

k  Solution 


fm 
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Reference  to  Tables  6  and 
7  show  that  LAOMO  =  -33530  „  . 

+  22.93T,  LMOAO  =  59910  a'\  ••.f«r.cu« 

-  27.11T  and  C  =  740S0  p.s  /V"  c.c»rund,„ 

Joules  respectively!  These  V//K\  r.  .< 

numerical  parameters  pro-  2mt  •  /7a  TV"  ‘L*  i52S,  / 

vide  a  means  of  defining  /II  '  / 

the  free  energy  of  the  ///  f  \  \\ 

compound  (1/7)  (At203-Mg0)  /  sv  \  \\  /  S'J 

(SP) .  The  free  energy  of  sr  ///  $  \  \\  /  7 

formation  of  this  compound  V/y  1  \  \  / / 

from  the  stable  forms  of  sf.c.s  \  \  \ c‘%f\sr >T"' 

alumina  and  magnesia  \#7  \  7j  \ 

(e.g.,  (1/S)  At20,(C)  and  c  *  fl  '•  \ 

(l/2)MgO(P)  is  readily  v  .  «  H  v  .■«.  * ..  »...  ^  _  ... 

obtained  by  adding  the  A,( 

quantities  0.714AOAOSC  and 
0.286  MOMOSP  to  the  ex¬ 
pression  for  GSP  in  Equa¬ 
tion  (11).  The  result,  itoor 

-5280  -  0.460T  Joule/g.at. 
is  then  multiplied  by 
seven  and  displayed  in 
Table  8  as  the  calculated 
free  energy  change  in 
Joules  for  the  formation 
of  Mg0-Ai203(SP)  from  c 

A^2  0j(C)  and  MgO(P)  i.e.,  A‘o.* «.»  i  . . 

-36961  -  3. 22T  Joules.  Figure  X9.  Calculated  Isothercjl  Sections  in  the  !!j0  S30  5-Fe0  4,9O5  57i 

It  is  also  compared  with  -Mo.«°o.6  sy*,e“- 

the  observed  value  (20) 
of  -35560  -  2.05T  Joules. 

In  summary.  Tables 
5-8  and  Figures  4-15 

illustrate  how  the  descriptions  of  the  elemental  edge  binary  systems  are  em¬ 
ployed  to  synthesize  the  quasi-binary  joins  of  interest  within  ternary  systems 
which  are  composed  of  the  elemental  edge  binaries. 

The  stability  of  the  computed  "quasi-binary  joins"  must  be  considered 
within  the  frame  work  of  the  approximations  noted  in  the  earlier  discussions  of 
"FeO"  (WO)  and  Fe203(F0).  In  particular,  each  of  the  computed  "joins"  shown  in 
Figures  4-15  can  be  undercut  by  more  stable  interactior.s  along  intersecting  joins 
in  the  same  ternary.  Such  interactions  are  not  likely  across  the  joins  between 
CO,  AO,  MO  and  HO.  However  joins  in  which  FO  and  KO  are  partners  may  be  sus¬ 
ceptible  to  such  intersections  by  competing  pairs  of  compounds  or  solution 
phases  within  the  elemental  ternary  system  which  lie  on  intersecting  joins! 

Such  interactions  must  be  considered  in  computing  the  most  stable  configuration 

The  extension  of  the  computations  into  quasi-ternary  systems  is  ill¬ 
ustrated  in  Figures  16-21  for  selected  isothermal  s  ctions  of  the  H0-C0-A0, 
MO-HO-AO  and  MO-CO-AO  systems  over  a  range  of  temperatures.  These  computations 
are  carried  out  in  a  straightforward  manner  emploj ing  the  methods  and  defini¬ 
tions  detailed  earlier.  Development  of  such  sections  as  a  means  for  guiding 
experimental  establishment  of  isothermal  sections  in  such  systems  is  probably 
the  most  useful  application  of  the  present  framework.  Indeed  very  little 
experimental  data  exist  with  which  to  compare  the  predicted  results.  Only  the 
isothermal  section  of  Greshkovich  and  Subican  (30)  in  the  MO-CO-AO  system  at 
1973K  which  is  shown  in  Figure  21  could  be  fountT’for  comparison  with  the  com¬ 
puted  results.  The  current  framework  can  also  be  employed  to  predict  interest¬ 
ing  compositions  for  eutectic  solidification  or  precipitation  reactions  on 
aging  after  high  temperature  solution  treatment.  The  periclase  solid  solution 
in  the  MO-HO-AO  and  MO-CO-AL  systems  displayed  in  Figures  18-21  exhibits  an 
extensive  range  of  stability  at  high  temperature  which  is  sharply  restricted 
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At  lower  temperatures. 

In  both  systems  the  com¬ 
plex  spinel  phase  limits 
the  range  of  stability  of 
the  periclase  structure. 
Under  such  conditions, 
solution  treatment  at 
high  temperatures 
(i.e.,  2S00K)  followed 
by  low  temperature  aging 
(i.e.,  below  1300K)  could 
be  employed  to  strengthen 
the  periclase  matrix. 

8.  Summar 


A  data  base  has 
been  presented  for 
describing  the  thermo¬ 
chemical  properties 
and  the  phase  diagrams 
for  a  group  of  six 
oxides.  Treatment 
of  these  oxides  as 
components  of  quasi¬ 
binary  joins  in  true 
ternary  systems  has 
been  employed  to 
develop  a  quasi -binary 
"p^roximation  which  can 
then  be  used  tc  compute 
higher  order  systems 
for  purposes  of  phase 
diagram  computation 
and  prediction  of 
.hermochemical  pro¬ 
perties. 


sswm 
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Figure  20.  Calculated  Isothermal  Sections  in  the  MfQSOQ>s*Cr0  4 0Q  6* 
^10.4°0.6  System. 
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APPENDIX  I 

Quasibinary  Approximation  to  the  Free  Energy 
of  Solution  Phases  (Mass  Basis,  One 
Gram  Atom  or  One  Mole  of  Atoms)  (See  Figure  1) 

GL«(l-x-y)°G^+x0G^y°G£+RT[(l-x-y)«.n(l-x-y)+x*nx+yi!.ny]+EGLtx,y,T]  (A-l) 

EGL[x,y ,T]*x(l-x-y) (1-y)'1 [ (l-x-y)LIJ*xLJI] +y (1-x-y) (1-x) "^ [ (1-x-y) LIK+yLKI] 

+xy(x+y)  *1[xLJK+yLKJ]  +xy( 1-x-y) TRNL  (A.2-) 

where  LIJ,  LJI,  LIK,  LKI,  LJK,  LkJ  and  TRNL  are  all  temperature  dependent  func¬ 
tions.  If  one  considers  the  vertical  plane  which  cuts  GL  and  contains  the 
compositions  x»xi,  y»o  and  x=x:n,  yal -r„,  with  a  quasi  binary  definition  of  the 
excess  free  energy  as 

EGL-5(l-5)l(l-5)L1+5L2]  (A- 3) 

where  ?-(x-xJ)  (x^-xj)'1,  y»5(l-x*)  and  (l-x-y)“(l-x£) (1-5) 

The  explicit  relationship  between  L,  and  L2  and  the  terms  LIJ,  LJI,  LIK,  LKI, 
LJK,  LKJ  and  TRNL  shown  in  Equation  (A- 2) can  be  obtained  by  equating  the  ternary 
and  quasi  binary  descriptions  shown  in  Figure  1  at  5-0,  1/3,  2/3  and  1.  The 
result  is  given  by  the  following  equations: 

L1*2M-N,  L2-2N-M,  P-(x*-xJ),  Q»(l-xi),  S-(l-x*)  (A-4) 

M^O.667  (l-0.333S)‘1(xi+0.333P)Q[Q0.667LIJ+(xi+P0.333LJI)]  (A- 5) 

Nj»0.333(l-0.667S)  '1(xi+0.667P)Q[Q0.333LIJ+(xl+P0.667LJI)]  (A-6) 

M2“0. 222 (Q-0. 333F) *1QS [Q0 . 667LIK+S0 . 333LKI)  (A- 7) 

N2«0.222(Q-0.667)‘1QS[Q0.333LIK+S0.667LKI]  (A-8) 

Mj*0 . 333 (xJ+0 . 333Q) "1S(xi+0 . 333P) { (xJ+0 . 333P)LJK+S0 . 333LKJ]  (A-9) 

N3«0.667(xi+0.667Q) '1S(xJ+0.667P) [(xi+0.667P)LJK+S0.667LKJ]  (A-10) 

M4-0.222TRNL  QS(xi+0.333P)  (A-ll) 

N4-0.222TRNL  QS(xJ+0 . 667P)  (A-12) 

Ms-RT[0.667QJ.n0.667Q+(xi+0.333P)in(xi+0.333P)+0.33Stn0.333S)  (A-13) 

Ns-RT[0.333Qin0.333Q+(xJ+0.067P)!.n(xi+0.667P)+0.667SS,n0.667S]  (A-14) 

M6-0.667Qx];(QLIJ+x];LJI)  (A-l£) 


-23- 


QUASIBINARY  AND  QUASITERNARY  OXIDE  SYSTEMS  -  I 


N6«0.333QxJ(QLIJ+xjLJI)  (A_lfi) 

M?-0 . 333SX* (SLKJ+x*LJK) 

(A-17) 

N.-0.667Sx*(SLKJ+x*LJK)  ..  10, 

MgH).667RT(x**nx*+Q4nQ)+0.333RT(x*Anx*+S4ixS-1.910)  (A.19) 

N8«0.333RT(xiEnxi+QEnQ-1.910)+0.667RT(x)lEnx)k+S£nS)  (A-20) 

M«4.S(M,.M2+M3+M4+Ms-M6-M7-M8)  (A-21) 

N-4.3(N1+N2+N3+N4+Ns-N6-N7-N8)  (A-22) 
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Larry  Kaufman  * 

ManLabs,  Inc.,  21  Erie  Street 
Cambridge,  Massachusetts  02139,  USA 


ABSTRACT.  A  data  base  is  being  developed  for  calculation  of  quasi-binary  and 
quasi-ternary  phase  diagrams  of  ceramic  systems.  The  initial  seg¬ 
ment  of  this  base  covers  combinations  of  Cr2Oj,  MgO,  Ai203,  Fe203 
Fe30^  and  "FeO",  Lattice  Stability,  Solution  and  Compound  Phase 
Parameters  were  derived  covering  the  liquid,  spinel,  periclase  and 
corundum  phases  in  systems  composed  of  these  oxides.  These  para¬ 
meters  permit  calculation  of  the  thermochemical  properties  and  phase 
diagrams  of  combinations  of  these  oxides  in  good  agreement  with 
limited  experimental  data.  This  base  has  been  extended  to  include 
Si02  and  CaO  by  carrying  out  an  analysis  of  the  thirteen  binary 
systems  composed  of  the  above  mentioned  systems  with  silica  and 
calcium  oxide.  The  current  base  of  eight  oxides  can  thus  be  em¬ 
ployed  to  compute  fifty  six  ternary  systems.  Isothermal  sections 
in  the  MgO-Fe jO„-Cr2Oi  are  presented  between  2800K  and  1600K  illus¬ 
trating  the  usefulness  of  the  data  base. 

1.  Introduction 

The  utility  of  computer  based  methods  for  coupling  phase  diagrams  and 
thermochemical  data  for  metallic  systems  has  been  demonstrated  in  many  papers 
published  in  this  journal.  A  considerable  effort  is  being  applied  toward 
development  of  a  suitable  data  base  for  metallic  systems.  Numerous  examples 
have  been  published  showing  how  such  data  on  metallic  systems  can  be  synthe¬ 
sized  to  successfully  predict  the  thermochemical  properties  and  phase  relations 
in  ternary  metallic  systems.  Such  predictions  can  be  applied  practically  to 
explore  conditions  for  in-situ  eutectic  growth,  metallic  glass  formation, 
spinodal  decomposition,  metastable  equilibrium,  and  diffusionless  transforma¬ 
tions.  In  addition,  recent  development  of  SIALON  composites,  combining  silicon 
and  aluminum  nitrides  with  oxides  of  silicon,  aluminum,  magnesium  and  other 
metals,  has  provided  additional  motivation  for  developing  methods  for  predicting 
multicomponent  phase  diagrams  of  ceramic  systems.  Accordingly  present  plans 
include  further  expansion  of  the  current  data  base  to  cover  SijNt,  and  ALN, 

In  addition  to  the  present  work  (1)  other  efforts  described  in  CALPHAD  are 
currently  in  progress  in  order  to  develop  predictive  methods  which  can  be  used 
to  compute  multicomponent  phase  diagrams  for  ceramic  systems  [2,3). 


•This  work  has  been  sponsored  by  the  Air  Force  Office  of  Scientific  Research 
Bolling  AFB,  Washington,  D.C.,  under  Contract  F44620-76-C-0060. 
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2.  Description  Of  The  Thermochemical  System 
Employed  To  Describe  Solution  And  Compound  Phases 


The  method  utilized  for  describing  solution  and  compound  phases -is 
essentially  the  same  as  that  employed  earlier  (1)  except  that  some  symbolic 
usage  has  evolved  which  facilitates  data  handling  as  indicated  below.  The 
free  energy,  GL,  of  a  liquid  (solution)  phase,  L,  in  the  binary  system  I-J  is 
given  by  Equation  (1)  where  T  is  in  Kelvins,  x  is  the  atomic  fraction  of  J  and 
the  mass  basis  is  a  mole  of  atoms  (i.e.,  a  gram-atom). 

Gl  ■  (1-x)°g£  +  x°Gj  +  RT(x£nx  +  (l-x)*n(l-x))+  x(l-x)  [LIJ(l-x)  +  xLJI]  (1) 

In  this  expression,  “G^  and  °G^  are  the  free  energies  of  a  gram  atom 
of  pure  liquid  I  and  J,  R  =  8.314  J/g.  at  .°K,  and  LIJ  and  LJI  are  functions 
of  temperature  (but  not  of  composition).  The  composition,  x,  is  the  atom  frac¬ 
tion  of  J.  Similar  equations  can  be  written  for  each  solution  phase  with  "G^ 
replaced  by  °G|,  °g£,  etc.,  for  the  P*Periclase,  C=Corundum  phases,  etc., 
and  with  LIJ  replaced  by  PIJ  and  CIJ  etc.  In  the  present  symbolism,  one  uses 
the  symbol  DODOLP  to  designate  the  free  energy  of  pure  liquid  CaO  minus  the 
free  energy  of  pure  Periclase  form  of  CaO.  Reference  to  Tables  1  and  2  show 
that  DODOLP=39748-13. 723T  Joules,  and  that  the  liquid  phase  of  the  CaO-MgO 
system  (i.e.  DO-MO  with  i“D0,  j=MO  and  x^atom  fraction  of  MO)  is  described  by 
LDOMO“-20920  and  LMODO=-20920.  The  thermochemical  data  for  stable  phases  listed 
in  Table  1  was  taken  from  the  ManLabs-NPL  Data  Bank  (4).  The  free  energy,  G*, 
of  a  compound  phase  *s  defined  by  Equation  (2)  as, 

GT  -  (1-x*)  °G®  +  x*  °G®  ♦  x* (1-x*) [ (l-Xft)LIJ  ♦  x*  LJI-C [T] ]  (2) 

where  6  refers  to  the  base  phase  selected  for  the  compound  and  C [T ]  (which  is  a 
temperature  dependent  function)  is  the  compound  parameter  .  Table  3  lists  the 
compound  parameters  and  base  phases  assigned  to  compounds  in  the  systems  of 
interest.  Thus  the  free  energy  of  one  gram  atom  of  the  compound  Ca0*Cr203 
which  is  defined  as  (1/7)  (Ca0’Cr203)  or  DOo.* ,  6C0i>.7 1  *  is  defined  as 

g'1'-  0.286  °Gp0  *■  0.714  °G^0  ♦  (0. 286) (0. 714) [ (0 . 286) LDOCO+O . 714LCODO-C]  (3) 

Reference  to  Table  2  and  3  shows  that  LDOCO*LCODO=0  and  the  compound  parameter 
C  in  this  case  is  equal  to  101671.  Thus  the  free  energy  of  formation  of 
DOo.jnsCOo.7i'.  from  the  spinel  forms  of  DO  and  CO  is  equal  to  -20741  Joules.  In 
order  to  compute  the  free  energy  of  formation  of  DOo.* esCO0.7 1 «  from  the  stable 
P(Periclase)  form  of  DO  (i.e.  CaO)  and  the  stable  C (Corundum) form  of  CO 
(i.e.  Cr203)  it  is  necessary  to  include  the  lattice  stability  terms  DODOSP  and 
C0C0SC.  The  former  is  contained  in  Table  1  (DODOSP=DODOLP-DODOLS*28137-4.93T) 
while  the  latter  was  presented  earlier  (1)  (i.e.  COCOSC=11756-0.92T) .  With 
this  additional  information  the  free  energy  of  formation  of  DOo.z < sCO0.7 1 4 
from  the  P  form  of  DO  and  the  C  form  of  CO  can  be  computed  as  is  shown  in  Table 
4.  Here  the  free  energy  of  formation  of  (1/7) (Ca0-Cr203)  from  (l/7)CaO 
(i.e.  0.286DO)  and  (l/7)(Cr203)  (i.e.  0.714  CO)  with  DO  in  the  Periclase  form 
and  CO  in  the  Corundum  form.  Reference  to  Tables  4  and  S  discloses  that  there 
are  differences  between  the  calculated  \alues  and  the  assessed  "Data  Bank 
Values"  which  are  largely  based  on  direct  thermochemical  measurements.  However, 
it  should  be  pointed  out  that  there  are  several  instances  of  "internal 
inconsistencies"  in  the  "Data  Bank  Values".  Thus  for  example  the  "assessed 
values"  in  Table  4  suggest  that  Ca0-Fe203  decomposes  into  Fe203  and  2Ca0.Fe203 
and  that  3Ca0*At203  decomposes  into  CaO  and  Ca0-A£203  at  300K.  Moreover 
the  "assessed  values"  in  Table  5  suggest  decomposition  of  mullite 
(3A£203 • 2Si02 )  into  Si02  and  A£203  at  low  temperatures.  Since  none  of  the 
above  reactions  are  actually  observed  it  must  be  assumed  that  these 
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TABLE  1 

SUMMARY  OF  LATTICE  STABILITY  PARAMETERS 
(All  units  in  Joules  per  gram  atom  (mole  of  atoms) , T in  Kelvins) 

P  «  Periclase,  C  =  Corundum,  S  =  Spinel,  X  =  Crystobalite , 

T  *  Tridymite,  H  =  hexagonal  (a  quartz),  R  =  trigonal  (B  quarts) 

AO  «  (1/S)A1203,  CO  =  (l/S)Cr203,  MO  »  (l/2)MgO,  WO  *  (l/2)FeO 
FO  -  (l/S)Fe203,  HO  =  (l/7)Fe304, 


SOSOLX*=  (1/ 3 )Si02 (Liquid)  -  (l/3)Si02  (Crystobalite) 
SOSOLT  «  (l/3)Si02 (Liquid)  -  (1/3) Si02  (Tridymite) 

SOSOTX  «  (l/3)Si02 (Tridymite)  -  (l/3)SiO£  (Crystobalite) 
SOSOLX  =  SOSOLT  +  SOSOTX 

SOSOLX  =  334V  -  1.674T 

SOSOLT  =  3933  -  2.008T 

SOSOTH  =  858  -  0.795T 

SOSOHR  =  251  -  0.293T 

SOSOLC  =  SOSOLP  =  SOSOLS  =  -2.092T 

DODOLP  =  39748  -  13.723T 

DODOLC  =  14644  -  10.878T 

DODOLS  «  11611  -  8.786T 

DODOLX  -  -  1.674T 

DODOLT  *  -  2.008T 


AOAOLX  ■  COCOLX  =  MOMOLX  =  WOWOLX  =  HOHOLX  =  FOFOLX  =  -1.674T 
AOAOLT  *  COCOLT  =  MOMOLT  =  WOWOLT  =  HOHOLT  =  FOFOLT  *=  -2.008T 

*  These  differences  specify  the  free  energy  of  one  phase 
(i.e.  liquid)  minus  the  free  energy  of  the  second  phase 
(i.e.  Crystobalite)  of  a  given  compound.  • 
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TABLE  2 

Quasibinary  Solution  Parameters  for  Oxide  Systems 
(All  units  in  Joules  per  gram  atom  (mole  of  atoms),  T  in  Kelvins) 


DO-AO  (l/2Ca0-l/SAl203) 

LDOAO  =  -2259S 
CDOAO  =  102926 
PDOAO  =  102926 
LAODO  =  -2259S 
CAODO  =  102926 
PAODO  =  102926 
DO-CO  (l/2Ca0-l/SCr„0,) 

nuflco  =  0  1  * 

CDOCO  =  83680 
PDOCO  =  83680 
LCODO  =  0 
CCODO  =  83680 
PCODO  *  83680 
DO-MO  (l/2CaO-l/2MgO) 
EBOMO  =  23012 
PDOMO  =  41840 
LMODO  =  29288 
PMOPO  =  68199 
DO-FO  (l/2Ca0-l/2Fe,0,) 
EBUFO  =  -20920  2  3 

CDOFO  =  104600 
PDOFO  =  104600 
LFODO  «=  -20920 
CFODO  =  104600 
PFODO  =  104600 
DO -HO  (l/2Ca0-l/7Fe,0.) 
LDOHO  =0  34 

SDOHO  =  1171S2 
POOHO  =  1171S2 
LHODO  =  0 
SHODO  =  117152 
PHODO  ■  117152 
DO -WO  (l/2CaO-l/2FeO) 
LDOWO  =  -7113 
PDOWO  *  30962 
LWODO  =  -20501 
PWODO  =  17573 


S0-D0  (l/3SiO,-l/2CaO) 
0<xdo<_0.40 

LSODO  “  0 

XSODO  =  83680 

TSODO  =  83680 

LDOSO  =  -393631+104. 6T 

XDOSO  =  -309951+104. 6T 

TDOSO  =  -309951+104. 6T 

0-40<xdq<1.0 

LSODO  =  -195476+41. 8T 
PSODO  =  -111796  +  41 .8T 
LDOSO  =  -100416+41. 8T 
PDOSO  =  -16736+41. 8T 
SO-AO  (l/3Si02-l/5Al,0,) 
TSMO  =  38493  2  3" 

LSOAO  =  122173 
XSOAO  =  122173 
TSOAO  =  122173 
LAOSO  =  28451 
CAOSO  =  1 1 2 1 3i 
XAOSO  =  112131 
TAOSO  =  112131 
SO- CO  (l/3Sift  -l/SCr,0T) 
“LSotO  =  681'9i>  2  3 

CSOCO  *  151879 
XSOCO  *  151879 
LCOSO  *  68199' 

CCOSO  =  151879 
XCOSO  =  151879 
SO -MO  (l/3SiO?-l/2MgO) 

LSOMO  =  -24267+24. 27T- 
XSOMO  -  59413*  24. 27T 
TSOMO  =  59413+24. 27T 
LMOSO  =  -229283+68. 62T 
XMOSO  =  -145603+68. 62T 
TMOSO  =  -145603+68. 62T 
0.40<xMO<l. 0 

LSOMO  *  -106274+42. 01T 
PSOMO  =  19246+42. 01T 
LMOSO  =  -106274+42. 01T 
PMOSO  -  19246+42. 01T 


SO-FO  (l/3Si02-l/SFe203) 

LSOFO  =  S4392 
CSOFO  =  138072 
XSOFO  =  138072 
TSOFO  =  138072 
LFOSO  =  163176-83. 68T 
CFOSO  =  246856-83. 68T 
XFOSO  *  246856-83. 68T 
TFOSO  =  246856-83. 68T 
SO -HO  (l/3SiO, - l/7Fe.O . ) 
OKX'HO<0.70  2  3  4 

LSOHO  =  80542-10. 46T 
XSOHO  =  164222-10. 46T 
TSOHO  =  164222-10. 46T 
LHOSO  =  213384-125. 52T 
XHOSO  =  297064-125. 52T 
THOSO  =  297064-125. 52T 

'0-70£xho^1-0 

LSOHO  =  0 
SSOHO  =  83680 
LHOSO  =  247902-130. OT 
SSOHO  =  331582-130. OT 
SO -WO  (l/3SiO,-l/2FeO) 

o<xwo<o.so 

LSOWO  =■  29288 
XSOWO  =  112968 
TSOWO  =  112968 
LWOSO  =  324678-225. 94T 
XWOSO  =  408358-225. 94T 
TWOSO  =  408358-225. 94T 
0.50<XWO<1.0 

LSOWO  =  162339-112. 97T 
PSOWO  =  246019-112. 97T 
LWOSO  =  191627-112. 97T 
PWOSO  =  275307-112. 97T 
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TABLE  3 

SUMMARY  OF  COMPOUND  PARAMETERS 

(All  units  in  Joules  per  gram  atom  (mole  of  atoms),  Tin  Kelvins) 


Compound 

Stoichiometry/ 

Base 

Compound  Parameter 
(Joules/g.atom) 

l/9(2Ca0-Fe203) 

DO0.44S 

FO0-S»S 

P 

9623  -  3.43T 

l/7(Ca0-Fe203) 

DO0.286 

F00.7Ti4 

S 

46861  -  1S.48T 

l/32(Ca0-6Al203) 

D00.063 

A0Q„<g37 

C 

16736  +  13.18T 

1/12 (CaO- 2A1203) 

DO0.167 

A00.S33 

C 

15481  +  10.46T 

l/7(Ca0-Al203) 

Dt)0.286 

A00-7ni4 

s 

38493 

l/ll(3Ca0-Al203) 

L00.546 

A00..«54 

p 

16318  +  13.51T 

l/7(Ca0-Cr203) 

D00.286 

C00-“il4 

s 

101671 

l/5(CaO-SiOz) 

D00.40 

S00_©0 

p 

-62383  +  49.29T 

l/7(2Ca0-Si02) 

DO0.571 

SO0-4t29 

p 

-56902  +  43.93T 

l/21(2Si02-3Al203) 

SO0.286 

AO0-7rH4 

c 

41840  +  4.73T 

l/S(MgO-SiOz) 

M00.40 

S00.®0 

p 

-66735  +  47.28T 

l/7(2Mg0-Si02) 

MO0.S71 

S00.4L29 

p 

-64434  ♦  41.13T 

l/7(Si02-2Fe0) 

S00.429 

W0O.S.71 

p 

217986  -  120. ST 

MO  «  (l/2)MgO,  DO  - 

(l/2)CaO 

,  FO  * 

(l/5)Fe203, 

AO  *  (1/5)A1203 

SO  *=  (l/3)Si02,  WO  * 

(l/2)FeO,  C  = 

Corundum,  P 

*  Perida  se 

S  *  Spinel 


inconsistencies  represent  cumulative  errors  which  can  occur  in  assessment  of 
thermochemical  data  when  the  competition  between  adjoining  phases  (i.e.  the 
phase  diagram)  is  not  considered,  iln  the  present  process  such  competition  is 
considered  as  a  natural  part  of  the  evaluation  method. 

3.  Discussion  off  ^Computed  Phase  Diagrams 

Binary  phase  diagrams  between  CaO  and  Si02  with  the  oxides  Ai;C3, 
Crj03  MgO,  FeO,  Fe20j  and  FejOi,  are  provided  in  the  compilation  series 
prepared  by  Levin,  Robbins  and  McMtedie (5).  Although  some  of  the  phase  diagrams 
are  sketchy  and  conflicting  experimental  data  are  presented  tht  results  have 
been  used  in  concert  with  the  experimental  thermochemical  data(4)  and  prior 
lattice  stability  values(l)  to  select  the  values  for  th*-  solution  and  compound 
parameters  which  are  shown  in  Tables  2  and  3.  Figures  1-13  are  the  binary 
phase  diagrams  which  can  be  generated  from  the  results  displayed  in  Tables  1-3. 
The  calculated  Atj03-Ca0  binary  diagram  in  Figure  1  is  in  relatively  good 
agreement  with  the  experimental  diagrams  shown  in  Figures  231,232,2295,2296  and 
4308  of  reference  S.  Moreover,  Ta&lle  4  shows  that  the  computed  free  energies 
of  formation  of  compounds  in  this  binary  system  are  in  relatively  good  agree¬ 
ment  with  the  "Data  Bank  Values"  derived  from  thermochemical  data.  The  solu¬ 
tion  phases  in  this  system  are  approximated  by  rpgular  solution  descriptions 
(i.e.  LDOAO=LAODO  etc.).  Figure  2  sShows  the  computed  Cr20}-Ca0  phase  diagram 


L.  Kaufman 


TABLE  4 


SUMMARY  OF  CALCULATED  AND  ASSESSED  (4)  FREE  ENERGIES 

OF  COMPOUND  FORMATION  (JOULES,  °K) 

DATA  BANK  AG  VALUES 

(2/9)CaO  + 

(l/9)Fe203-  l/9(2CaO-Fe203) 

-5117 

at 

300K 

AG  - 

-2714-0. 313T  (calculated) 

-7627 

at 

1400K 

(1/7)C*0  + 

(l/7)Fe203-*  l/7(Ca0-Fe203) 

-3112 

at 

3,0  OK 

AG  * 

-3808  +  0.977T  (calculated) 

-4284 

at 

1400K 

(l/32)CaO 

+  (6/32)Al203-  (l/32)(Ca0-6Al203) 

-- 

-- 

AG  «= 

-739-0. 957T  (calculated) 

-- 

-- 

(l/12)CaO 

+  (2/12)Al203-  (1/12) (Ca0-2A1203) 

-1916 

at 

300K 

AG  = 

-1101-1. 928T  (calculated) 

-5381 

at 

1600K 

(l/7)CaO  + 

(1/7)A1203-  (l/7)(Ca0«Al203) 

-3096 

at 

300K 

AG  - 

-1248-2. 219T  (calculated) 

-6720 

at 

1600K 

(3/ll)CaO 

+  (1/11)A1203-  (l/ll)(3Ca0-Al203) 

*1632 

at 

300K 

AG  « 

-947-2. 15ST  (calculated) 

-5786 

at 

1600K 

(l/7)CaO  + 

(l/7)Cr203-  (1/7) (Ca0-Cr203) 

-- 

-- 

AG  - 

-4300-2. 067T  (calculated) 

-- 

which  compares  well  with  Figure  37  of  reference  5.  Table  2  discloses  that 
this  system  is  computed  with  regular  solution  parameters  and  an  ideal  liquid 
phase  (i.e.  ID0C0=LC0D0=0) .  Figure  3  displays  the  computed  CaO-MgO  system 
which  is  obtained  with  assymetric  excess  free  energies  (LDOMO/LMODO)  but  ideal 
mixing  entropies  (dL/dT=0) .  The  latter  agrees  with  Figure  229  of  reference  5. 
Figures  4-6  contain  the  calculated  Ca0-Fe203,  CaO-FejO*  and  CaO-FeO  phase 
diagrams  which  are  in  good  agreement  with  those  shown  in  Figures  43-48  of 
reference(S) .  Reference  to  Table  2  shows  that  the  CaO-FejOj  and  CaO-FejOi. 
cases  are  approximated  by  regular  models  (i.e.  the  latter  as  an  ideal  liquid) 
while  the  CaO-FeO  case  is  assymetric  (LDOWO/LWODO)  with  ideal  entropies  of 
mixing  for  solution  phases.  It  should  be  noted  that  the  computed  Ca0-Fe203 
diagram  shown  in  Figure  4  is  metastable  since  Fe203  decomposes  into  FesO,  and 
O2  above  1760K  under  atmospheric  conditions. 

In  contrast  to  the  results  encountered  in  analysis  of  the  CaO-based 
systems,  reference  to  Table  2  shows  that  analysis  of  the  Si02  -based  binary 
systems  proved  to  be  more  difficult.  In  this  group,  only  the  ■Bi02-Cr203  binary 
calculated  in  Figure  9  could  be  approximated  by  a  regular  solution  model. 

The  resulting  phase  diagram  compares  well  with  the  experimental  diagram  shown 
in  Figure  332  of  reference  5.  Figure  P  shows  the  computed  Si02-Aft20j  diagram 
which  is  in  good  agreement  with  the  exi- :rimental  results  which  are  given  in 
Figures  313,4374  and  4375  of-Reference  5.  This  system  is  described  by  assy¬ 
metric  solutions  (LSOAO/LAOSO)  and  ideal  solution  entropies  (dL/dT  etc.=*0). 
Moreover,  examination  of  Table  5  shows  that  the  computed  free  energy  of  mullite 
(3Al203*2Si02)  from  the  stable  forms  of  Ai20i  and  Si02  is  in  relatively  good 
agreement  with  the  values  listed  in  the  Data  Bank(4).  The.  latter  are  deter¬ 
mined  by  assessment  of  thermochemical  data.  The  computed  Si02-Fe20}  diagram 
shown  in  Figure  11  is  metastable  above  1760K  due  to  decomposition  of  Fe20s 
into  FesOi,  and  02  at  one  atmosphere.  This  diagram  is  in  accordance  with 
Figure  80  of  reference  5  and  is  computed  with  assymetric  so.  tion  models 
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TABLE  5 


SUMMARY  OF  CALCULATED  AND  ASSESSED  (4)  FREE  ENERGIES 

OF  COMPOUND  FORMATION  (JOULES, °K) 

DATA  BANK 

AG 

VALUES 

(l/S)CaO  + 

(l/S)Si02-  (l/S)(CaO-Si02) 

AG  ■ 

-19791-2. 402T  (calculated  at. 

300K) 

-17970 

at 

300K 

AG  « 

-20808-1. S47T  (calculated  at 

1800K) 

-18146 

at 

1800K 

(2/7)CaO  + 

(l/7)Si02-  (1/7) (2CaO-Si02) 

AG  * 

-1848S-0.954T  (calculated  at 

300K) 

-19430 

at 

300K 

AG  « 

-19212-0. 343T  (calculated  at 

1800K) 

-19539 

at 

1800K 

(2/21)SiOz 

+  (3/21)Al20,-  (l/21)(2Si02 

•3A1203) 

AG  « 

-706-1. 2S4T  (calculated  at  300K) 

+397  at  300K 

AG  - 

-1191-0. 846T  (calculated  at 

1900K) 

-2013 

at 

1900K 

(l/5)MgO  ♦ 

(l/S)Si02-  (l/S)(MgO-Si02) 

AG  * 

-6464-1. 864T  (calculated  at 

300K) 

-7273 

at 

300K 

AG  - 

-7481-1. 009T  (calculated  at 

1800K) 

-5799 

at 

1800K 

(2/7)MgO  ♦ 

(l/7)Si02-  (l/7)(2MgO-Si02) 

AG  - 

-8088-0. 216T  (calculated  at 

300K) 

-8301 

at 

300K 

N 

<3 

-8815  +  0.39ST  (calculated  at  1800K) 

-7761 

at 

1800K 

(l/7)Si02 

+  (2/7)Fe0^  (l/7)(SiO2-2Fe0) 

AG  * 

-7373  ♦  1.413T  (calculated  at  300K) 

-6845 

at 

300K 

AG  * 

-7849  +  1.880T  (calculated  at  1400K) 

-4067 

at 

1400K 

(i.e.  LSOFO/LFOSO)  with  finite  excess  entropies  of  solution(i.e.  dL/dT/O)  as 
shown  in  Table  2. 

Analysis  of  the  remaining  binary  systems  i.e.  Si02-Ca0,  SiQ2-Mg0, 
Si02-Fe30,.  and  Si02-Fe0  could  only  be  carried  out  by  segmenting  the  description 
of  the  liquid  phase  as  is  illustrated  in  Table  2.  Thus,  in  the  last  of  the 
above  named  cases,  the  excess  free  energy  of  the  liquid  phase  (i.e.  the  last 
term  in  Equation  ,1)  is  set  equal  to  x(l-x)  [(l-x)29288+x(324678-225.94T)J  where 
O^x^O.S  and  x=xwo=atom  fraction  of  WO.  In  the  range  0.5<x<1.0  the  last  term 
is  set  equal  to  x(l-x) ((l-x) (162339-112. 97T)+x(191627)-112.97T) j.  Note  that 
both  expressions  yield  the  identical  result,  44246-28. 24T  at  xwo=xson0.5. 

The  calculated  Si02*Fe0  phase  diagram  shown  in  Figure  13  is  in  relatively  good 
agreement  wi*h  the  observed  diagram  in  Figure  80  of  reference  5  except  that 
liquid  monotectic  is  reported(S)  at  1970K  and  45  w/o  FeO.  Examination  of  Table 
5  shows  that  the  computed  free  energy  of  formation  of  Fayalite  (2Fe0*Si02)  from 
the  stable  forms  of  Si02  and  FeO  is  in  relatively  good  agreement  with  the 
"assessed  values"(4)  over  the  temperature’ range  of  its  stability.  In  the  Si02- 
FejOi,  case  the  thermochcmical  description  of  the  liquid  phase  shown  in  Table  2 
is  divided  into  two  segments  corresponding  to  the  composition  range  0<xho£0.70 
(i.e.  H0=(l/7)Fc30<,)  and  a  second  covering  the  range  0.70<xjiq£1.0.  TKe  computed 
binary  phase  diagram  shown  in  Figure  12  compares  favorably  with  that  shown  in 
Figure  87  of  reference  5  except  that  the  observed  liquid  monotictic  is  reported 
(5)  to  be  at  1950K  and  90  w/o  FejOi,. 
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w/o  CaO 

Af203  13  26  37  48  58  67  76  85  93  CaO 


w/o  CaO  — t-* 

Cr2°3  9  18  27  37  47  57  67  78  89  CaO 
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w/o  CaO 


Figure  7.  Calculated  SiQ  33jOfl  667-CaQ  sOQ  5  Phase  Diagram 


w/o  At^Oj 


S0  MULLITE«(l/21)(2SiO2-2AE2O3)«SO0+286AO0>714  A0 
Si0/.333°0.667  At0.4°0.6 

Figure  8.  Calculated  SiQ#  jjjOg^-AR^O^  Phase  Diagram 
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v/o  Fe203  — - 

SiO,  IS  29  41  52  62  71  79  86  94  Fe2°3 


'  1  1  1  1 
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_  H(aQuartz)  :+  P 

S*  5  0 

R(BQuartz)  f  P(Periclase) 


Sl0.333°0.667 


Fe0.4°0.6 


Figure  11.  Calculated  SiQ  333O0  667-FeQ  4Oq  6  Phase  Diagram 

(Note:  FeQ  40Q  &  Decomposes  into  FeQ  429°o  571  and  1/2  °2  above  1760K) 

w/o  Fe.O. _ _ 


iv  1  •  1  I 


Ll+L2  \ 

_ _ iiu _ 

X(Crystobalite)  +  L 


T(Tridymitej  +  S(Spinel) 


H(aQuartz)  +  S 


R(BQuartz)  +  S(Spinel) 


0.3330 . 667 


T+L  L(Liquid) 
.  )  S+L 


,w0. 429^0. 571 


Figure  12.  Calculated  Si0>333O0>667-Fe0>429O0#571  Phase 
Diagram 
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l 
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S°0.429WO0.S7r(1/7^SiO2*2FeO^  W° 

Fe0  5®0 . 5^ 


0.333v0.667 

Figure  13.  Calculated  sio.33300.667*Fo0.50O.S  phasc  Dia8ram 


In  the  remaining  Si02-Ca0  and  Si02-Mg0  cases.  Table  2  shows  that  the 
liquid  thermochemistry  is  described  in  two  composition  ranges,  from  1.0>XgQ>p.6 
and  from  0.6>xso^.°-  The  resulting  computed  binary  diagrams  shown  in  Figures 
7  and  10  compare  favorably  with  those  shown  in  Figures  237,  2331,  2303  ?n^  766, 
2331  of  reference  S.  In  these  cases  the  greatest  discrepancies  are  encountered 
again  in  the  position  of  the  liquid  monotectic  where  the  experimental  locations 
are  at  1970K  and  27  w/o  CaO  in  the  Si02-Ca0  case  and  at  1960K  and  31  w/o  MgO  in 
the  SiOi-MgO  case.  Comparison  of  the  computed  and  "assessed"  free  energies  of 
compound  formation  for  SiOj-CaO  and  SiO* -MgO  compounds  listed  in  Table  5  shows 
that  in  both  cases  excellent  agreement  is  observed  for  the  (1/7) (2Mg0-5i02 ) 
(i.e.  MOo.5  7iSOo.4  2-.)  and  (1/7) (2Ca0-Si02) (i.e.  DOo.s 7 1  SO0.4 2 9)  cases.  However, 
while  the  computed  low  temperature  results  for  the  (1/5) (MgO-Si02) 

(i.e.  MO0.4SO0.5)  and  (1/S) (Ca0-Si02 )  (i.e.  DO0.4SO0.6)  agree  with  the  "assessed" 
values,  the  computed  high  temperature  values  are  in  poor  agreement  with 
the  "assessed"  high  temperature  values.  Thus,  it  appears  that  while  the 
present  model  can  be  readily  applied  to  obtain  a  good  description  of  the  CaO 
based  systems  difficulity  is  encountered  in  its  application  to  some  of  the  Si02 
based  systems. 


4.  Calculation  of  Quasi-Ternary  Phase  Diagrams 

The  objective  of  developing  the  present  data  base  for  characterizing 
quasi-binary  ceramic  systems  is  to  provede  a  means  for  computing  stable  and 
meta  stable  equilibria  in  quasi-ternary  and  higher  order  ceramic  systems. 
Several  examples  of  such  systems  were  presented  earlier (1).  In  these  cases  the 
quasi-ternary  systems  were  computed  by  synthesis  of  the  edge  binary  character¬ 
istic  for  solutions  on  the  basis  of  the  Koehler  equation.  TJ:”s  for  the 
(l/2)MgO- (l/7)Fe304- (l/5)Cr203  (i.e.  MO-HO-CO)  case  the  free  energy  of  the 
liquid  is  defined  as 
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GL-(l-x-y)  G^q  +  xG|j0  +  *  RT((l-x-y)in(l-x-y)+x*nx+yiny) 


♦  x(l-x-y) (1-y)  [(l-x-y)LMOHO+xLHOMO] 


+y(l-x-y)(l-x) 


-1 [ (l-x-y)LMOCO+yLCOMO] 
+xy (x+y) ~ 1 [xLHOCO  +  yLCOHO] 
♦TRNL(l-x-y)xy 


(4) 


where  x  is  the  atom  fraction  of  HO,  y  is  the  atom  fraction  CO  and  TRNL  is  a 
ternary  interaction  parameter.  Similar  equations  can  be  written  for  the 
Periclase(P) ,  Corundum(C),  Spinel (S)  and  other  solution  phases  by  replacing 
the  L’s  and  TRNL  by  P's  and  TRNP  etc.  If  TRNL=0  then  all  of  the  terms  in 
Equation(4)  have  been  explicitly  defined(l).  Thus  in  this  approximation  the 
free  energy  of  ternary  solution  phases  are  synthesized  from  their  binary  edge 
components  by  using  Kohler's  Equation  with  a  zero  ternary  interaction  parameter. 
The  free  energy  of  ternary  compound  phases  is  defined  along  specific  joins 
As  an  example,  we  consider  the  spinel  compound  phase  in  the  system  MO  HO-CO 
(see  Figure  14).  This  compound  runs  between  the  Fe3Oi,'MgO  and  CrjOj'MgO, 
compositions  designated  (Table  7(1))  as  (Feo.j#s00.»o»)  •  (Mgo.i  <•  jOo.i  » 3)  and 
Cr«.j »(0(i 4 2 1 )  •  (Mgo  j  „  jOo.j  *  j) :  If  x  is  the  fraction  MO  and  y  is  the  fraction  HO 
then  the  free  energy  of  the  spinel  (SP)  compound  CO,,  l,MO  l-HOri  ,MO  . 
given  by  Equation(S)  as  u'x*'  **  U*x»)  **  is 

GSP  -  (l-x-y)Gpn  ♦  xGyn  ♦  YG*  ♦  RT[yEny+(l-x-y)in(l-x-y)  -  (l-x)in(l-x)  ] 


JCO  '  A“MO  '  7 “HO 
♦  [1-  + 


where 


AGg  -  0,AGa  -  (l-xi)xj[(l-xj)  LCOMO+xJlmOCO  -  CA] 


AG 


B  -  (l-x*)x*l(l-x*)  LHOMO  ♦  x*  LMOHO  -  Cfi] 


(5) 


(6) 


and 


x«xj  +  yp  ;  p  =>  (x*  -  xj)/(l-x,),xi 


0  and  C. 


38910+16. 74T 


0.286,  x*  *  0.286,  p  ■  0 
(From  Table  7  of  reference  1) 


(7) 


Since  all  of  the  parameters  required  to  calculate  the  free  energy  of 
the  liquid,  spinel,  periclase  and  corundum  solution  phases  in  the  M0-H0-C0 
system  as  well  as  .the  free  energy  of  the  spinel  compound  phase,  the  entire 
ternary  phase  diagram  can  be  computed.  The  results  are  shown  as  a  series  of 
isothermal  sections  between  2800K  and  1600K  in  Figures  14  and  15.  These 
results  can  he  employed  to  locate  melting  points  and  liquid  decomposition 
products  thus  providing  a  guide  for  structural  studies.  Moreover  these  sections 
can  be  employed  to  locate  interesting  solid  state  decomposition  reactions  which 
can  be  employed  to  produce  strengthening  reactions  attainable  via  heat  treat¬ 
ment.  One  such  example  in  the  present  computations  is  the  decomposition  of 
the  periclase  field  indicated  on  cooling  from  1850K  to  1600K  in  Figure  15. 

Thus  an  "alloy"  composed  of  O.7MO-0.25HO-0.U5CO  (i.e.  88.4  mol  percent  MgO-9.1 
mol  percent  FejO»-2.5  mol  percent  Crj03)  would  be  expected  to  be  strengthened 
by  solution  treatment  at  18S0K  and  subsequent  aging  at  1600K. 


; 

i 


L.  Kaufman 
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CALCULATION  OF  QUASIBINARY  AND  QUASITERNARY 
OXYNITRIDE  SYSTEMS  -  III* 
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ManLabs,  Inc.,  21  Erie  Street 
Cambridge,  Massachusetts  02139  USA 

(This  paper  was  presented  at  CALPHAD  VIII,  Stockholm,  Sweden, May  1979) 


ABSTRACT.  A  data  base  is  being  developed  for  calculation  of  quasi-binary  and  quasi-ternary 
phase  diagrams  of  ceramic  systems.  Initially  the  base  was  restricted  to  the 
following  oxides:  CrjOj,  MgC,  Fe^j,  Fe^°4»  Fe0>  A12°3’  Ca0  an<*  S*°2"  Lattice 

Stability,  Solution  and  Compound  Phase  Parameters  were  derived  covering  the 
liquid,  spinel,  corundum,  periclase,  crystobalite,  tridymite  and  quartz  phases 
which  appear  in  the  binary  systems  composed  of  these  oxides.  These  parameters 
were  selected  in  accordance  with  the  observed  thermochemical  properties  and 
phase  diagrams  of  these  binary  systems.  The  usefulness  of  this  analytical 
description  was  illustrated  by  computing  isothermal  sections  for  ternary 
systems  composed  of  these  oxides.  The  above  mentioned  data  base  has  now  been 
expanded  to  cover  combinations  of  the  oxides  with  silicon  nitride,  Si^N4,  and 
aluminum  nitride,  AIN  thus  permitting  computation  of  oxynitride 
SIALON  ceramic  systems  which  are  of  current  practical  interest.  The  current 
analysis  provides  isothermal  sections  in  the  Si  N^-AljO  -AIN  and  SiC^-A^Oj- 

SijN^  systmes  as  well  as  computed  values  for  the  free  energy  of  formation  of 

the  beta  sialons  in  the  Si^-Al^-AIN  system  and  the  X  phase  in  the  Si02-Al203- 

SijNj  system.  The  former  results  are  compared  with  independent  calculations 

reported  by  Doerner  et.al.  at  CALPHAD  VIII  and  published  in  the  current  issue 
of  the  CALPHAD  Journal. 


1.  Introduction 

The  utility  of  computer  based  methods  for  coupling  phase  diagrams  and  thermo¬ 
chemical  data  for  metallic  systems  has  been  well  documented  in  many  papers  published  in 
this  journal.  A  considerable  effort  is  being  applied  toward  developing  an  extensive  data 
base  for  metallic  systems.  Recently  similar  efforts  have  begun  with  the  aim  of  developing 
similar  data  bases  for  ceramic  systems  (1-4).  Recent  development  of  SIALON  composites, 
combining  silicon  and  aluminum  nitrides  with  oxides  of  silicon,  aluminum,  magnesium  and 
other  metals  has  provided  additional  motivation  for  developing  methods  for  predicting 
multicomponent  phase  diagrams  of  ceramic  systems.  The  present  paper  describes  work  con¬ 
ducted  at  ManLabs  during  the  past  year  and  compares  the  results  obtained  independently  at 
the  Max  Plank  Institute  in  Stuttgart  (4). 

2.  Description  of  the  Thermochemical  System 
Employed  to  Describe  Solution  and  Compound  Phases 

The  method  utilized  for  describing  solution  and  compound  phases  is  the  same  as 
that  employed  earlier  (1,2)  employing  some  symbolic  usage  which  facilitates  data  handling 
as  indicated  below.  The  free  energy,  GL,  of  a  liquid  (solution)  phase,  L,  in  the  binary 
system  I-J  is  given  by  Equation  (1)  where  T  is  in  Kelvins,  x  is  the  atomic  fraction  of  J 
and  the  mass  basis  is  a  mole  of  atoms  (i.e.,  a  gram-atom). 

Gl  «  (l-x)*Gj  ♦  x°Gj  +  RT(x*nx  ♦  (l-x)£n(l-x))+  x(l-x) [LIJ(l-x)  ♦  xDI]  (1) 

*  This  work  has  been  sponsored  by  the  Air  Force  Scientific  Research,  Bolling  AFB, 
Washington,  D.C.,  under  Contract  F44620-76-C-0060. 
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In  this  expression,  °Gb  and  °Gb  are  the  free  energies  of  a  gram  atom  of  pure 

liquid  I  and  J,  R  =  8.314  J/g.at.  °K,  and  LIJ  and  UI  are  functions  of  temperature  (but 
not  of  composition).  The  composition,  x,  is  the  atom  fraction  of  J.  Similar  equations 

can  be  written  for  each  solution  phase  with  8Gb  replaced  by  °G?,  8g9,  etc.,  for  the 

P=Periclase,  C=Corundur  phases,  etc.,  and  with  LIJ  replaced  by  PIJ  and  CIJ  etc.  In  the 
present  symbolism,  onr  uses  the  symbol  SNSNLB  to  designate  the  free  energy  of  pure  liquid 
SijN^  minus  the  free  energy  of  the  beta  form  of  Si^.  Reference  to  Tables  1  and  2  show 

SNSNLB=33949- 12 . S IT  Joules/g.at. ,  and  that  the  liquid  phase  of  the  SN-AO,  (l/7)Si3N4-(l/5)Al203 

system  (i.e.  SN-AO  with  i=SN,  j=AO  and  x=atom  fraction  of  AO)  is  described  by  LSNA0=29288 

and  LA0SN=29288.  The  free  energy,  G  ,  of  a  compound  phase  i,  J  is  defined  by  Eq.  (2)  as 

■  (l-x#)  x. 


G*  =  (1-xJ  °G®  ♦  x,  #G®  ♦  x#(l-x„)[(l-xJLIJ  +  xt  UI-C[T]]  (2) 

where  0  refers  to  the  base  phase  selected  for  the  compound  and  C[T)  (which  is  a  temperature 
dependent  function)  is  the  compound  parameter.  Table  2  lists  the  compound  parameters 
and  base  phases  assigned  to  compounds  in  the  systmes  of  interest.  Thus  the  free  energy  of 
one  gram  atom  of  the  compound  AlN'A^Oj  which  is  defined  as  (1/7)  (AlN’A^Oj)  or 

^0.286^0.714  is  defined  as 


Ct=  0.286  "G^  ♦  0.714  8G®Q  ♦  (0.286) (0.714) [ (0.286)LANAO+0. 714  LAOAN-C]  (3) 


Reference  to  Table  2  shows  that  LANAO=LAOAN=29288  and  the  compound  parameter  C  in  this  case 
is  equal  to  146440.  Thus  the  free  energy  of  formation  of  ANQ  286AO0  714  from  the  beta 

SijN4  forms  of  AN  and  AO  is  equal  to  -23922  Joules.  In  order  to  compute  the  free  energy 
of  formation  of  A1Q  286A0Q  714  ^rom  tbe  stable  N(hexagonal)  form  of  AN  (i.e.  AIN)  and  the 
stable  C(Corundum)  form  of  A0  (i.e.  A^Oj)  it  is  necessary  to  include  the  lattice  stability 

terms  ANANBN  and  A0A0BC.  The  former  is  contained  in  Table  1  (ANANBN=ANANLN-ANANLB=42593- 
0.80T)  while  the  latter  AOAOBC=AOAOLC-AOAOLB=23640+2.30T.  Kith  this  additional  information 
the  free  energy  of  formation  of  ANQ  286*%  714  ^rom  the  N  fonn  of  ***  and  t*ie  C  ^orm  A0 

can  be  computed  as  is  shown  in  Table  3.  Here  the  free  energy  of  formation  of  (1/7) 
(A1N*A1203)  from  (1/7)A1N  (i.e.  0.286AN)  and  (1/7)  (AI^)  (i.e.  0.714  AO)  with  AN  in  the 

hexagonal,  N,  form  and  AO  in  the  Corundum  form  as  illustrated  by  Equation  4. 

AG  =  0.286  Cg]^  -  'G^)  ♦  0.7? '  (8gJq-  8gJq)-23922  (4) 

«  0.286  (ANANBN)  ♦  0.714  (A0A03C)  -  23922  (5) 


*  5138  ♦  1.418T  Joules/g.at.  (6) 

Tables  1  and  2  summarize  the  lattice  stabilities,  solution  phase  and  compound 
phase  parameters  for  the  quasi-binary  oxynitride  systems  of  interest  while  Table  3  displays 
the  free  energy  of  formation  of  the  compound  phases  which  form  in  the  subject  binary  systems. 
The  quasi-binary  systems  computed  on  the  basis  of  Tables  1  and  2  in  combination  with  the 
lattice  stability  values  presented  earlier  are  shown  in  Figures  1-7  which  are  all  computed 
at  one  atmosphere  pressure!  Thus  these  diagrams  represent  metastable  equilibria  at  tempera¬ 
tures  above  the  decomposition  points  of  AN  and  SN.  This  matter  will  be  taken  up  in  Section 
3.  Reference  to  Table  3  shows  that  several  metastable  phases  are  described  in  the  AN-AO 
system.  These  phases  designated, ANQ  706AOQ  2g4(G)  and  ANQ  286AOo  714^  display  positive 

free  energies  of  formation  of  1.121T  Joules/g.at.  and  5138+1. 418T  Joules/g.at.  respectively. 
Indeed  the  latter  was  discussed  in  Equations  (3-6)  above.  In  equation  2  it  was  shown  that 
the  free  energy  of  formation  of  this  ANQ  2g6AOQ  714(B)  phase  from  AN(B)  and  A0(B)  was 

negative  i.e.  -23922  Joules/g.at.  However  Equations  (3-6)  which  compute  the  free  energy 
of  formation  of  this  phase  from  AN(N)  and  A0(C)  yield  the  positive  value  in  Table  6.  Thus 
this  phase  does  not  appear  on  the  AN-AO  diagram  shown  in  Figure  3.  Reference  to  Figure  3 
shows  that  neither  the  above  mentioned  AN  706A0  2Q4(G)  nor  the  ^0  286A®0  714 ^  phase 


-46- 


TABLE  1 


SUM4ARY  OF  LATTICE  STABILITY  PARAMETERS 
(All  units  in  Joules  per  gram  atom  (mole  of  atoms),  T  in  Kelvins) 


AO  *  (1/S)A1203,.S0  =  (l/3)Si02>  MO  =  (l/2)MgO,  AN  =  (1/2)A1N,  SN  =  (l/7)Si3N4 
P  *  Periclase,  C  =  Corundum,  S  =  Spinel,  X  =  Crystobalite,  T  =  Tridymite 
H  *  a  quartz,  R  =  8  quartz,  B  =  beta  SijN^,  N  =  hexagonal  AIN,  L  =  Liquid 

SNSNLB*  =  (1/7)  Si3N4  (Liquid)  -  (1/7)  Si^fbeta) 


ANAjNLN 

— 

:  (1/2)  AIN  (Liquid)  -  (1/2)  AIN 

(hex)  etc. 

SNSNLB 

= 

33949 

-  12. SIT 

ANANLN  =  42S93 

- 

13.31T 

SNSNLN 

= 

-  13.31T 

ANANLB  = 

- 

12. SIT 

SNSNLP 

= 

-  12.S5T 

ANANLP  = 

- 

10.21T 

SNSNLC 

= 

-  10.21T 

ANANLC  =  . 

- 

10.21T 

SNSNLX 

= 

-  1.67T 

. 

AN  AN  I.X  = 

- 

1.67T 

SNSNLS 

= 

-  8.37T 

ANANLS  = 

- 

8.37T 

MOMOLP 

= 

47698 

-  1S.40T 

(reference  1) 

M0M0LB 

= 

-  12. SIT 

A0A0LB  = 

- 

12. SIT 

M0M0LN 

= 

-  13.31T 

AOAOLN  = 

- 

13.31T 

SOSOLC 

= 

-  2.09T 

(reference  2) 

AOAOLC  =  23640 

- 

10.21T  (reference  1) 

SOSOLB 

= 

-2.09T 

SOSOLN 

= 

-2.09T 

SOSO LX 

= 

3347 

-1.674T 

(reference  2) 

*  These  differences  specify  the  free  energy  of  one  phase  (i.e.  liquid)  minus  the  free 
energy  of  the  second  phase  (i.e.  beta)  for  a  given  compound. 

is  stable.  Indeed  these  phases  are  the  unstable  counterphases  for  the  stable 

SNfl  412^0  588^  31,(1  the  SN0  980^0  020^  phases  which  appear  in  the  SN-AN  system  shown 

in  Figure  2.  This  definition  of  the  metastable  ANQ  706AOQ  2g4(G)  31,(1  Mq  286A00  714^B^ 

phases  are  required  for  calculation  of  the  SN-AO-AN  ternary  system  described  in  Section  5. 
Tables  1  and  2  permit  computation  of  the  melting  points  of  these  metastable  compounds.  In 
particular  subtraction  of  equation  2  from  equation  1  provides  the  description  of  the  free 
energy  of  fusion  of  compound  phase  ?  when  x=x#. 

Thus 

AG™  -  GL-Gf*(l-x4)(°G!;  #G’]*x,[’cJ-#G®]+RT(x.enx,»(i-x,)in(l-x,))*x.(l-x#)C  (7) 

The.  free  energy  change  on  fusion  for  the  ANg  286A%  714^  phase  with  6=B  and  C=146440  is 
given  by 

GL-GB  *  0 . 286 (ANANLB) +0.714 (A0A0LB) -4 . 977T+29904  (8) 

«  29904-17. 487T  Joules/g.at  (9) 

yielding  a  melting  point  of  1710*K  which  lies  below  the  stable  2200*K  eutectic  in  Figure  3. 
A  similar  calculation  can  be  performed  for  the  ANg  7Q^AO(  2g4(G)  phase  with  8=N  and 

062760  from  Table  2. 
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TABU-  2 


SUMMARY  OF  SOLUTION  AND  COMPOUND  PARAMETERS 
(All  units  in  Joules  per  gram  atom  (mole  of  atoms}.,  T  in  Kelvins) 


Solution  Parameters 


LANSO  =  LS0AN  =  LSNAN  =  LANSN 

=  LANAO  = 

LAOAN  =  LANMO  = 

LMOAN  =  29288 

LSNSO  =  LS0SN  =  LSNAO  =  LAGSN 

=  tMOSN  = 

LSNMO  =  29288 

NANSO  =  NS0AN  =  XANS0  =  XS0AN 

=  TANSO  r 

TSOAN  =  HANSO  = 

HSOAN  =  RAN SO 

= 

P.S0AN  =  125520 

BSNAN  =  BANSN  =  N'SN'AN  =  NAN3N 

=  125520 

NANAO  =  NA0AN  =  CANAO  =  CA0AN 

=  125520 

NANMO  =  NM0AN  =  PANMO  =  PM0AN 

=  125520 

a 

BSN50  =  BS0SN  =  XS0SN  =  XSN30 

=  TSNSO  = 

:  TSOSN  =  HANSO  = 

HSOAN  =  RANSO 

= 

RSOAN  =  125520 

BSNAO  =  BA0SN  =  CSNAO  =  CAOSN 

=  125520 

BSNMO  =  BMOSN  =  PSNMO  =  PMOSN 

=  125520 

Compound  Parameters 

Compound 

Parameter 

Compound 

Name 

Stoichiometry 

Stability  Base  Joules/g.at. 

(1/50) (7Si3N4-(l/2)AlN) 

B 

^0.980^0.020 

Stable 

B 

83680  ♦  4.184T 

(1/17) (Si^-SAIN) 

G 

^0. 412^0. 588 

Stable 

N 

96232 

(1/17) (6A1N  •  A1203) 

G 

AN0.706A°0.294 

Metastable 

N 

62760 

(1/  7) (AIN-  A1203) 

B 

^0.286^0.714 

Metastable 

B 

146440 

(1/17) (AIN  5A1203) 

SP 

^0.118^0.882 

Stable 

S 

120081  -  12.678T 

(1/12)(5A1N-Mg0) 

Q 

"‘o.sss^o.ie? 

Stable 

N 

104600—  8.368T 

(l/50)(7Si3N4-(l/3)Si02) 

b 

SN0.980S°0.020 

Stable 

B 

71128  ♦  4.184T 

(l/10)(Si3N,-Si02) 

R 

^O^OO^O.SOO 

Stable 

B 

115060*  25.104T 

(1/50)  (7Si3N4- (1/5)A1203) 

B 

SN0.980A00.020 

Stable 

B 

6S270  ♦  4.184T 

(l/SO)((l/2)MgO-7Si3N4) 

B 

M00.020SN0.980 

Stable 

5 

112968  ♦  4.184T 

Here 

GL-GG  =  G.706(ANANLN)*0.294(AOAOLN)-S.036T+I3027  (10) 

*  43098-18. 346T  Joules/g.at.  (11) 

yielding  a  melting  point  of  2349*K  which  is  unstable  since  it  lies  within  the  two  phase 
t+N  field  shown  in  Figure  3. 

Finally  mention  should  be  made  of  the  SnQ  gggAN0  020(B),  Sn0  980^°0  020 ^ 

S»O.MOA00.0!0"»  “■*  "°0.  020SnO  980^®^  compound  phases  shown  in  Tables  2  and  3  and  in 

Figures  2,5,6  and  7.  These  "compound"  phases  are  defined  for  mathematical  convenience 
only  so  that  the  extension  of  the  Si.N^  or  SN(B)  phase  into  a  ternary  system  can  be  made 

along  the  lines  used  previously  for  metallic  systems (5).  Examination  of  the  free  energy 
of  formation  of  these  compounds  show  very  small  negative  values  which  could  be  representative 
of  a  dilute  solution  phase. 
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TABLE  3 


SttWARY  OF  CALCULATED  FREE  ENERGIES  OF  COMPOUND  FORMATION 
(JOULES,  ®K) 


REACTION 

(7/SO)Si3N4(B)  ♦  U/100)AlN(N)-(l/S0)(7Si3N4-O/2)AlN)(B) 
(l/17)Si3N4(B)  ♦  (S/17)A1N(N)  ->(1/17)  (Si^- SAIN)  (G) 
(6/17)AlN(N)  +  (l/l'7)Al203(C)  -(1/17) (6A1N-A1203) (G) 
(1/7)A1N(N)  «■  (1/7)A1203(C)  ->(1/7)  (AlN'AljOj)  (B) 

(1/17)A1N(N)  +  (3/17)Al203(C)  ->(1/17)  (AIN- 3A1203)  (SP) 
(S/12)A1N(N)  ♦  (l/12)MgO(P)  -(1/12) (SAlN-MgO) (Q) 

(7/S0)Si3N4(B)  ♦  (l/lS0)SiO2(R)-(l/S0)(7Si3N4-(l/3)SiO2)(B) 
(l/10)Si3N4(B)  ♦  (l/10)Si02(R)-(l/10)(Si3N4-Si02)(R) 
(7/50)Si3Ni4(B)  >  (l/S)Al2O3(C)-(l/S0)(7Si3N4-(l/S)Al2O3)(B) 
(l/100)MgO(P)  >  (7/SO)Si3N4(B)-(l/SO)((l/2)MgO-7Si3N4)(B) 


AG 

-213  -  0.098T 
-2234  ♦  0.326T 
-  1.121T 
S138  ♦  1.418T 
-SIS  -  0.2S9T 
-2510  ♦  0.81ST 
-720  -  0.102T 
16498  >  4. 97 IT 
-234  -  0.036T 
-686  -  0.140T 


3.  Calculation  of  the  Melting  and  Decomposition  Temperatures  of  AIN  and  hl3"4 


Examination  of  Figures  1-7  shows  that  the  decomposition  temperatures  for  AIN 
and  Si3N4  are  computed  to  be  2850K  and  2129K  respectively  at  one  atmosphere .Ihese  tempera¬ 
tures  are  below  the  computed  melting  points  of  AN  and  SN  which  are  defined  by  ANANLN  and 
SNSNLB  in  Table  1  at  3200K  and  2714K  respectively.  Calculation  of  the  melting  and 
decomposition  temperatures  of  AN  and  SN  was  carried  out  along  the  lines  developed  previously 
(1,6-12)  by  defining  Equations  1  and  2  for  the  liquid  and  A1Q  50NQ  50  phase  in  the  Al-N 

system  and  the  liquid  and  SiQ  42gNo  571  P^ase  *n  t*ie  Si'N  system.  Table  4  shows  the 

results  in  terms  of  Lattice  Stability  Parameters  (4.1),  Solution  Phase  Parameters  (4.2), 
Compound  Parameters  (4.3)  and  Free  Energy  changes  on  Fusion  4.4.  The  latter  provides  the 
basis  for  the  values  of  ANANLN  and  SNSNLB  in  Table  1.  The  equations  displayed  earlier  in 
connection  with  the  melting  point  of  ANq  286A®0  714^  ani*  ^0  706^*0  294^  aru*  derived  by 

subtraction  of  equations  1  and  2  (i.e.  Eqs.  7-11)  are  used  in  defining  the  results  in  (4.4). 
The  Table  4  description  can  also  be  employed  to  compute  the  free  energy  of  formation  of 
AN  and  SN  from  nitrogen  and  liquid  or  solid  forms  of  aluminum  and  silicon.  The  results 
shown  in  4.5  are  in  good  agreement  with  the  assessed  Data  Bank  values  (8).  These  results 
yield  the  predicted  one  atmosphere  decomposition  points.  Finally,  the  solubility  of 
nitrogen  in  liquid  silicon  and  liquid  aluminum  can  be  computed  (4.6)  and  are  found  to  be 
in  good  agreement  with  observed  values  (9,10), 

4.  Calculation  of  the  'i02-Al203-Si3N4  Ternary  System 

The  isothermal  section  at  2200K  in  the  Si02-Al203-Si3N4  system  was  computed  and 

is  shown  in  Figure  8.  The  ternary  liquid(L),  corrundum(C)  and  B(beta)  solution  phases  were 
described  by  means  of  Kohler's  equation  with  zero  values  for  the  ternary  interaction  para¬ 
meters.  The  mullitc(M)  and  silicon  oxynitride(R)  compounds  were  defined  as  shown  in 
Figure  8  with  the  counterphase  paraws-ters  CB=0  and  CAB=0  along  the  lines  previously  set 
forth  in  Equations  4-6  of  reference  12) .  In  addition  a  ternary  X  phase  was  included 
corresponding  to  SOQ  2QA0  65SfJ0  15  comPoun<*  parameter  for  this  phase 
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Figure  4.  Calculated  A1Q  5NQ  g  -MgQ  500>5  Phase  Diagram  at  one 
atmosphere 
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Figure  S.  Calculated  Mg-  ^0.  £  -Si.  42gNo  571  pfiase  DiBgram  at 
one  atmosphere 

(i.e.  C=47,700)  shown  in  Figure  8  is  defined  on  the  basis  of  Table  VIII  of  reference  12 
which  defines  the  free  energy  of  a  ternary  compound  by  means  of  Equation  (12)  as 

G*  *  0.20  ‘G^q+0.65  “G^+0.15  °GgN+(0.20)  (0.6S)  (0.85)'1  [0.20LSOAO+0.65LAOSO-C] 

♦  (0.20) (0.15) (0.35)”* [0.20LSOSN+0. 15LSNS0-C]* (0.65) (0. IS) (0.80)”*[0.65LADSN+0. 15LSNAO-C] 

(12) 

Since  LAOSO=28451  and  LSOAO=38492  from  reference  (2) 

G*  «  0.20  ‘G^q+0.65  °G^0+0.15  °GgN-9460  Joules/g.at.  (13) 

or  the  free  energy  change  for  the  reaction 

0.067Si02(x)+0.13Al203(C)+0.0214SijN4(BM0.067Si02-0.13Al203.0.0214Si3N4](X)  (14) 

AG  *  0.20(S0SOCX)+0. 15(SNSNCB)-9460  =  -3700-0. 261T  Joules/g.at.  (15) 


S.  Calculation  of  the  SijN^-Al^-AlN  Ternary  System 

The  ternary  isothermal  section  for  the  SN-A0-AN  system  shown  in  Figure  6  was 
calculated  along  lines  similar  to  those  described  in  Section  4.  The  solution  phases 
covered  included  the  liquid,  L,  corundum,  C,  hexagonal,  N,  and  Beta,  B,  phases  described 
by  Kohler's  equation  with  zero  valued  ternary  interaction  parameters.  In  addition,  the 

spinel  ^0.118^0. 882 ^SP^’  SN0. 412^0. S88_AN0.706A00. 294 ^  and  SN0.980A00.020” 

AN0  286AO0  714^  compound  phases  are  included.  The  spinel  compound  phase  is  described 
with  a  SNq  ijgA0Q  gs2 CSP)  counterphase  with  6=S,  CB  and  CAB  equal  to  zero.  The  base 

phases  and  the  compound  parameters  for  the  B  and  G  compound  phases  are  shown  in  Figure  9. 
The  prior  discussion  in  Section  2  described  the  metastable  counterphases  ANQ  7q^ao0  294 ^ 

and  ANq  286AOo  714  W*  These  metastable  phases  are  described  with  base  phases  of  N  and  B 

respectively  and  with  compound  parameters  which  are  62760  and  146440  Joules/g.at. 
respectively  as  shown  in  Table  2.  Panels  (a)  through  (c)  in  Figure  9  show  the  computed 
equilibrium  between  the  N/L,  B/L  and  G/L  pairs  at  2200K.  The  B/L  equilibrium  shown  in 
Panel  (b)  was  computed  with  a  value  of  CAB=-83862  in  order  to  extend  the  B  (beta  SijN4) 
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TABLE  4 


SIM1ARY  OF  ALUMINUM-NITROGEN  AND  SILICON-NITROGEN  SYSTEM 
ANALYSES  FOR  ESTIMATION  OF  THE  MELTING  TEMPERATURES  AND 
HEATS  OF  FUSION  OF  AIN  and  SijN  (Joules,  g.at.  °K) 


V=Vapor,  L=Liquid,  D=Diamond  Cubic,  A=fcc,  B=bcc 
» 

1.  Lattice  Stability  Parameters 


ALALLB=  628-  6.69T  (reference  6) 

ALALLA=10711-11.S1T  " 

ALALBA=10083-4.82T  " 


NNNNLB=  879-17S7T 

NNNNBV=-3682+53.64T 

NNNNLV=-2803+36.07T 


SISIBD=44350-19.54T  (reference  S) 

SISILD=50626-30.00T  " 

SISILB=  6276-10. 46T  " 


2.  Solution  Phase  Parameters 


LALNN=-S0208,  LNNAL=-912112+2S6.48T,  LSINN=14644,  L=-674879+23S.14T 

3.  Compound  Phase  Parameters 

(Al0  5Nq  5)  ;  C* 167360+ 18. 4 IT  base=bcc 

(Si0  429Nq  571)  ;  C=125520+31.38T  base=bcc 

4.  Free  Energy  Change  on  Fusion 


A10.SN0.S: 

AG  *  0. 5(628-6. 69T)+0. 5(879-17. 57T)-5.77T+0. 25(167360+18. 41T) 

*  42593-13. 31T  =  ANANLN (Table  1)  ;  Melting  Point=3200°K 

Si0.429N0.S71: 

AG  *  0. 429(6276-10. 46T)  +  0.571(879-17, 57T)  -  5. 68T+0. 245(125520+31. 38T) 

*  33949-12. SIT  =  SNSN LB (Table  1)  ;  Mel;ing  Point=2713°K 


5.  Free  Energy  of  Formation 
0.5  A1(L)  +  0.25N2(V)  +A10>5N0<5 

AG  =  0.50(r628+6.69T)  +  0.5(-3682+53.64T) 

♦  0. 25 (0. 50 (-50208) +0.50 (-912 1 12+256. 48T)- 167360-18. 41T) 

*  -164285  ♦  57.625T;  Decomposition  at  1  atmosphere  at  2850°K 

0.5  Al(fcc)  +  0.25N2(V)  -+  A1Q  SNQ  s 

AG  =  0.50(10083-4. 82T)  +  0.5(-3682+53.64T) 

♦  0. 25 (0. SO (-50208) +0.50 (-9121 12+256, 48T) -167360-18. 41T) 

*  -158929  +  51.87T 


TABLE  4  (concluded) 


0.429  Si(L)  *  0.285SN2(V)  -+  Si0<429N0>s71 

AG  -  0.429 (-6276+ 10. 46T)  +  0.571(-3682+S3.64T) 

♦  0.245(0.429(14644)  +  0.571(-674879+235.14T)-12SS20-31.38T) 

•  -128420+60. 32T;  Decomposition  at  1  atmosphere  at  2129®K 

0.429  Si(D)  +  0.2855N2(V)  -+  Si0429N0S71 

AG  *  0. 429(44350-19. S4T)  +  0.571 (-3682+53. 64T) 

+  0.245(0.429(14644)  +  0.571(-674879+23S.14T)-125520-31.38T) 

*  -106702+47. 44T 

6.  Calculation  of  Nitrogen  Solubility  in  Liquid  Aluminum  and  Silicon 
Al-N  (Solubility  of  nitrogen  in  liquid  aluminum) 

RTlnx  =  ALALBL  +  NNNNBL  -  LALNN 

+  0. 50 (0.5 (-50208) +0. 50 (-9121 12+256. 48T)- 167360- 18. 41T) 

*  -275559  +  79.18T  ;  x=0. 0001=0. Olatomic  percent  at  1773K 

Si-N  (Solubility  of  nitrogen  in  liquid  silicon) 

RTlnx  =  (0.429)  (0.S71)'1  SIS I BL+NNNNBL-LS INN 

+  0.429 (0.429 (14644)+0.S71 (-674879+235 . 14T) -125520-31 . 38T) 

*  -236701+69. 55T  ;  x=0. 00017=0. 017  atomic  percent  at  1673K. 


compound  phase  out  into  the  ternary  system  (3,4).  Panel  (d)  shows  a  composite  of  the  B/L, 
G/L,  N/L  and  S/L  equilibrium.  This  composite  constitutes  the  ternary  section  at  2200K  and 
should  provide  a  picture  of  the  phase  constitution  at  lower  temperatures.  Calculation 
of  ternary  sections  at  lower  temperatures,  (wher.  the  liquid  disappears)  can  be  readily 
carried  out  on  the  basis  of  the  parameter  values  developed  here.  It  should  be  noted  that 
the  isothermal  sections  shown  in  Figures  8  and  9  are  metastable  with  respect  to  the  one 
atmosphere  decomposition  of  Si^  at  2129°K  (see  Section  3) . 

6.  Calculation  of  the  Free  Energy  of  Formation  of  Beta  SIALONS  and  Comparison 
with  the  Computed  Results  of  the  Stuttgart  Group(4) 

The  foregoing  results  can  be  compared  directly  with  those  presented  earlier 
(3,4)  by  the  Stuttgart  group  with  respect  to  phase  equilibria  and  the  free  energy  of  forma¬ 
tion  of  the  beta  SIALONS.  Figures  b  and  9  should  be  compared  with  Figure  7  of  reference 
3  and  Figure  16  of  reference  4.  With  respect  to  the  former  comparison  Figure  8  covering 
the  SN-SO-AO  ternary  section  should  be  compared  with  the  upper  left  portion  of  the  SiC>2- 

AljOj-AlN-SijNj  square  in  Figure  7  of  reference  3  which  corresponds  to  the  SijNj-SiO-Al^ 

right  isosceles  triangle.  Here  the  X^  phase  shown  in  Figure  7  of  reference  3  corresponds  to 

the  current  X  phase  in  Figure  8.  Figure  9  should  be  compared  with  the  lower  right  isoceles 
triangle  Sij^-A^Oj-AlN  shown  in  Figure  7  of  reference  3.  The  present  diagram  includes 

the  "beta  SijNj"  phase  and  the  G  phase.  The  latter  corresponds  to  the  x7  phase  in 
reference  3.  The  phases  labeled  X2,X4>Xj  and  X^  in  reference  3  are  not  included  in  Figure 

9.  The  next  step  in  the  comparison  is  computation  of  the  free  energy  of  formation  of  beta 
sialons.'  This  calculation  is  detailed  in  Table  5  starting  with  the  expression  for  G® 
given  by  Equation  (16)  in  Table  5.  Reference  to  Equation  (16)  and  the  succeeding  relations 
through  Equation  (21)  shows  that  the  free  energy  of  the  beta  sialon  phase  which  projects 
into  the  Si3N4-Al203-AlN  ternary  from  the  SijN4  comer  toward  the  A1N'A1203  or  Al.jO.jN 

composition  can  be  defined  explicitly  in  terms  of  the  compound  parameter  for  the 


Calculated  Pairwise 


TABLE  S 


CALCULATION  OF  THE  FREE  ENERGY  OF  FORMATION  OF  BETA  SILICON 
NITRIDE  -  ALUMINUM  OXYNITRIDE  SOLID  SOLUTIONS 
(Joules,  °K) 


GB  *  z  °g|n  +  x  4GB0  +  y  ‘G^  ♦  (l-(y/(I-x)))AGA  +  (y/(l-x)  AGfi 

♦  RT(ylny  +  zlnz  -  (l-x)ln(l-x))  ♦  AG£  (from  reference  5)  (16) 

where  x=xAQ  and  y=yAN  and  zaz^l-x-y  as  in  Figure  9 


Since  the  solid  solution  runs  from  SnQ  g80AOQ  Q2Q  to  ANQ  2g6AOo  714  t*len  x*=0.714, 
xj=0.020  and  p=(x„-x|)/(l-x„)=2.426.  Thus 


x=x**yp  =  0.02  +  2.46y  (from  reference  5)  (17) 

z=l-x]  -  y(l+p)  =  0.98  -  3.426y  (18) 

with  AGa  =  (0 . 02) (0 . 98) (0 . 02LSNAO+0 . 98LAOSN -C] =0 . 0196 [29288-65270-4 . 184T] 

«  -705  -  0.082T  (19) 

and  AGb  =  (0.286) (0.714) [0.286LANAO+0.714LAOAN-C]=0. 2042(29288-146440] 

*  -23922  (20) 


since  the  compound  parameters,  C,  in  these  equations  refer  to  the  values  listed  in  Table  2 

for  SN0>980A00.020  3nd  ^0.286^0.714  resPectively 
Thus 

G®  •=  (0.98-3.426y)°G|N  ♦  (0.02+2.426y)#GB0  ♦  y  *G^ 

+  (l-y(0.98-2.426y)*1) (-70S-0.082T)  ♦  y(0.98-2.426y)"1(-23922) 


♦  RT(ylny+(0.98-3.426y) In (0. 98-3. 426y)-(0. 98-2. 426y) In (0. 98-2. 426y)) 

♦  CABy(0. 98-3. 426y) (0.98-2. 426y)-1  (21) 

The  free  energy  change  for  the  reaction 
(0.14-0.490y)Si3N4(B)+(0. 004+0. 485y)Al203(C)+0.5y  A1N(N)  •+ 

[(si3N4) (O.i4-0.409y) {A12°3} (0.004+0. * 85y) fA1N)0.5y^ (B)  (22) 


is 

AG  =  G®  -  (0.98-3.426y)°G|N  -  (0.02+2,426y)°GAO  -  y 
or  since  *GB0  -  *GAQ  =  AOAOBC=AOAOLC-AOAOLB=23640+2.30T 

from  Table  1  and  Reference  I  and 

*CAN  ’  =  ANANBN=42593-0. 80T  from  Table  1 


(23) 

(24) 

(25) 


-S8. 


TABLE  5  (concluded) 


Then  Equation  (23)  yields  the  following  result  when  CAB=-83862  J/g.at. 

AG  =  (0. 02+2. 426y) (23640+2. 30T)  ♦  y(42S93-0.80T) 

+  (l-y(0. 98-2. 426y)_1) (-705-0. 082T)  ♦  y(0.98-2.426y)-1 (-23922) 

♦  RT(ylny+(0. 98-3. 426y)ln(0. 98-3. 426y)-(0. 98-2. 426y) In (0. 98-2. 426y)) 

-83862y(0.98-3A26y)(0.98-2.426y)*1  (26) 

for  0<y<0.286.  Thus  Eq.  (26)  permits  explicit  calculation  of  the  free  energy  of  formation 
of  the  beta  silicon  nitride-aluminum  oxynitride  solid  solution  as  a  function  of  concentra¬ 
tion  along  the  join  between SN((l/7)Si3N4)  ;nd  ANQ  286A00  714^1/7  AlN-A^O^or 

(1/7)A1jNjO.  If  Y  is  the  mole  fraction  of  Al^N^O  along  the  Si^N^-Al^N^O  join  then  Y  and 

y  (the  atom  fraction  of  AN  in  the  SN-AO-AN  ternary  shown  in  Figure  9)are  simply  related  by 

Y  =  (0.98/0. 286)y+0. 02  =  3.426y+0.02  (27) 


metastable  ANQ  286A00  714  P^ase»  C=146440,  given  in  Table  2  and  the  value  for  CAB=-83862 

Joules/g.at.  which  was  used  to  compute  Figure  9.  Equations  21  through  26  show  how  the  free 
energy  of  formation  of  the  beta  sialon  phase  can  be  computed  as  a  function  of  composition 
when  the  compounds—  beta  Si^N^,  the  corundum  form  of  AljO^  and  the  hexagonal  form  of  AIN 

are  considered  to  be  reactants.  Thus  equations  (26)  and  (27)  can  be  employed  to  compute 
the  free  energy  of  formation  of  beta  sialons  according  to  equation  22  as  a  function  of 
composition.  Figure  10  shows  the  result  of  such  a  calculation  which  was  kindly  carried  out 
by  Dr.  H.L.  Lukas  (3,4)  as  a  function  of  the  mole  percent  AljOjN.  The  mole  percent  is  100 

times  the  mole  fraction,  Y,  given  in  Equation  27.  In  Figure  10,  the  curve  labeled 
L.  Kaufman  is  that  given  by  Equation  26  while  the  curve  labeled  P.  Doemer  et  al  is  the 
identical  free  energy  change  as  computed  by  the  Stuttgart  group  from  a  completely  different 
starting  point  using  an  entirely  different  model  for  the  liquid  and  beta  sialon  phase.  The 
differences  between  the  two  curves  are  less  than  4KJ/g.at.  (i.e.  1  kcal/g.at.)  which  is 
quite  remarkable.  If  the  values  of  CAB-33862  and  C=146400  Joules/g.at.  employed  in  Table  5 
to  develop  Equation  26  are  reconsidered, the  difference  between  the  two  curves  can  be 
reduced  to  less  than  1600  J/g.at.  (i.e.  400  cal/g.at.).  Thus  the  current  value  of  C=146440 
yields  a  positive  free  energy  of  formation  for  the  metastable  compound  ANQ  2g6AOQ  714 

as  noted  in  Section  2,  Equations  4-6  and  Tables  2  and  3.  Moreover  the  melting  point  of 
this  compound,  based  on  C=146,400  Joules/g.at.  is  1710K  according  to  Eq.  9.  If  the  value 
of  C  were  increased  to  C=167,360  Joules/g.at.  Equations  4-6  and  9  would  lead  a  free  energy 
of  formation  (from  hexagonal  AIN  and  A1203(C))  of  866+1. 418T  Joules/g.at.  as  compared  with 

the  current  value  of  S138+1.418T  listed  in  Table  3.  Thus  the  increase  in  C  would  make  the 
free  energy  of  AnQ  286A0Q  714  more  negative  than  it  is  currently,  but  still  positive. 

Finally  alteration  of  C  from  146,400  to  167,360  would  increase  the  melting  point  of 
AN0  286A00  714  frora  -•7i0°K  t0  1956°K»  which  is  still  below  the  2200K  eutectic  in  Figure  3. 

Thus  if  C=167,360  Joules/g.at.  for  ANQ  2g6AOQ  714  the  upper  curve  in  Figure  10  would  be 

sheared  downward  until  it  coincided  with  the  lower  curve  at  100  mole  percent  Al  0  N. 

*  O 

Conversely  one  could  inquire  as  to  the  accuracy  of  the  lower  curve  at  100  mole  percent 
AljOjN  since  the  latter  phase  is  unstable!  The  final  step  in  bringing  the  upper  curve 

close  to  the  lower  curve  is  decreasing  CAB  from  -83862  Joules/g.at.  to  -92048  Joules/g.at. 
The  effect  of  such  a  change  on  the  computed  phase  diagram  (along  with  the  aforementioned 
change  in  C  from  146400  to  167360)  is  shown  in  Figure  11.  The  net  result  is  an  increase 
in  the  field  of  stability  of  beta  SijN4, 
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T*  • 


6.  Conclusions 

The  current  extension  of  the  quasibinary 
synthesis  of  oxynitride  ceramic  systems  has  provided  a 
means  for  synthesizing  the  thermochemical  properties  and 
the  phase  diagrams  of  complex  multicomponent  systems.  The 
utilization  of  metastable  phase  descriptions  as  a  means 
for  assessing  phase  stability  has  been  shown  to  be  of 
value  in  deriving  information  on* the  stability  of  beta 
sialons.  The  current  description  has  also  been  applied 
to  compute  the  free  energy  of  formation  of  the  X  phase 
in  the  SiC^-SijN^-A^Oj  system  and  could  be  applied  to 

calculation  of  other  oxynitride  solid  solution  phases 
(i.e.  the  G  phase  in  the  Si^-AljOj-AlN  system).  It  i)> 

hoped  that  the  future  development  of  this  interesting 
and  inportant  CALPHAD  area  will  provide  future  oppor¬ 
tunities  to  perform  comparison  of  independent  comparisons 
along  the  lines  of  Figure  10  with  increasing  levels  of 
agreement ! 


a.M  N  00  <«  00  tO .09  <0.00  lOfl.io 

H01C-Z  AL303N 


Figure  10.  Comparison  of  Current 
Computed  Free  Energies  of  Forma¬ 
tion  of  Beta  Si.N,  Solid  Solu- 
3  4 

tions  at  2000K  with  those  due  to 
Doemer  et  al  (4). 


SN=(l/7)Si3N4 


f  AN=(1/2)A1N 


A0b(1/5)A1203 


Figure  11.  Recalculated  Composito  Equilibrium  in  the  (l/7)Si3N4-(l/5)Al20j-(l/2)AlN 
System  at  2200K  with  CAB=-92048  Joules/g.at.  and  C=167360  Joules/g.at. 
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UTILIZATION  OF  DATA  BASES  AND  COMPUTER  TECHNIQUES 
FOR  SOLVING  METALLURGICAL  PROBLEMS* 

by 
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ABSTRACT 

Documentation  of  specific  instances  in  which  thermochemical  data  have 
been  employed  successfully  to  economic  advantage  is  difficult  because  of  the 
secrecy  and  confidentiality  restrictions  which  of  necessity  limit  publication 
of  such  information.  Nevertheless,  such  examples  come  to  light  over  a  period 
of  years  and  have  been  properly  documented.  Moreover,  development  or  compu¬ 
tational  facility  can  be  assessed  as  providing  a  vehicle  for  evaluating, 
guiding  and  synthesizing  potential  economic  projects.  The  CALPHAD  group  is  an 
international  association  of  practicing  thermochemists  devoted  to  fostering 
technological  advances  in  thermochemical  knowledge  and  its  application  to  real 
problems.  The  fourth  meeting  of  the  CALPHAD  group  held  in  August  1975  at  the 
National  Bureau  of  Standards  in  Gaithersburg,  Maryland  provided  a  list  of 
examples  demonstrating  the  importance  of  phase  diagram  and  thermochemical  data 
in  developing  industrial  processes.  The  proceedings  of  this  meeting  were 
published  in  the  CALPHAD  Journal .  Succeeding  meetings  of  the  CALPHAD  group, 
as  well  as  subsequent  journal  papers  have  provided  many  nore situations.  The 
current  paper  will  review  these  cases  in  the  fields  of:  in  situ  columbium, 
nickel  and  cobalt  base  supcralloys,  sigma  phase  formation  in  high  temperature 
alloys,  multicomponent  semiconductors  for  electro-optical  devices,  SIALON 
ceramics, molten  salt,  sulfide  and  oxide  slags  for  ore  reduction  and  refining, 
high  temperature  battery  electrodes,  hardcnaMlity  calculations  in  alloy  steels, 
synthesis  of  hydrogen  storage  alloy  compounds,  grain  refinement  in  aluminum 
alloys,  metatectic  reactions  in  hafnium  alloys  and  the  growth  of  HEM  silicon 
solar  cell  material. 

Introduction 

.  The  symposia  and  publications  of  the  CALPHAD  group1'3  which  have 
taken  place  during  the  past  five  years  have  provided  many  examples  of  the 
application  of  thermochemical  data  bases  to  the  solution  of  industrial  problems 
The  CALPHAD  Journal,  published  regularly  by  Pergamon  Press  of  Oxford,  England 
provides  a  vehicle  for  timely  exposition  ■  of  information  and  instances  and 
such  applications.  Nevertheless,  it  is  useful  to  compile  additonal  illustrations 
from  time  to  time  in  order  to  demonstrate  the  wide  range  of  applications  and 
to  review  examples  presented  in  other  publications.  Accordingly,  the  present 
paper  will  cover  recent  cases  related  to  the  topics  listed  in  the  abstract 
which  demonstrate  how  data  bank  thermochemical  data  can  be  applied  to  the 
solution  of  real  problems. 


*  This  work  has  been  sponsored  by  the  Air  Force  Office  of  Scientific  Research 
Bolling  AFB,  Washington,  D.C.  under  Contract  F  44620-76-C-0060  and  the 
Metallurgy  Program,  Metallurgy  and  Materials  Section,  Division  of  Materials 
Research,  National  Science  Foundation,  Washington,  D.C.  under  Grant 
n\tP7fi-fl34A3  _ 
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1.  Calculation  of  c.uperalloy  Phase  Diagrams.  Today  there  is  strong  interest 
in  the  development  of  directionally  solidified  eutectic  alloys  for  the  blades 
and  vanes  of  the  gas  turbine  engine  which  are  being  subjected  to  increasingly 
higher  temperatures  and  stresses.  These  alloys  derive  their  superior  high 
temperature  strength  from  the  directional  solidification  of  an  alloy  of 
eutectic  composition  in  a  high  temperature  gradient  producing  aligned  high 
strength  microdimensional  fibers  or  lamellae  in  a  relatively  ductile  matrix. 
The  potential  for  increases  in  average  operating  metal  temperatures  of  over 
100#C  have  already  been  demonstrated  for  the  y/y'+S  and  carbide  reinforced 
superalloy  in  situ  composites.  Limitations  in  terms  of  oxidation  resistance, 
ductility,  and/or  castability  in  these  first  generation  alloys  indicate  the 
need  for  discovering  other  eutectic  alloy  systems.  These  may  be  ternary  or 
higher  order  systems,  so  there  is  'a  requirement  for  systematically  and 
efficiently  identifying  promising  multi-element  eutectics.  A  substantial 
number  of  publications  are  devoted  to  displaying  methods  for  computing  ternary 
isothermal  sections  of  superalloy  systems6'10.  Detailed  computations  of  the 
Cb-Al-Cr  system10  published  in  1973  and  shown  in  Figure  1  predicted  a  eutectic 
with  a  minimum  melting  point  at  1317°C  and  a  composition  of  41Cr-25Cb-34Al 
(wt.%)  within  the  region  of  the  Cr-Cb-CbAl3-Cr2Al3  composition  space.  This 
investigation  suggested  that  directional  solidification  of  an  alloy  of  this 
composition  should  result  in  an  in  situ  composite  having  a  CbAl3  reinforcing 
phase  in  a  matrix  of  64Cr-10Cb-26Al  (wt.%),  providing  a  combination  of  good 
high  temperature  strength  and  oxidation  resistance.  Subsequent  experimental 
investigation  of  this  system  by  M.K.  Thomas11  shown  in  Figure  2  verified  the 
eutectic  and  matrix  compositions  and  disclosed  a  meiting  point  at  ^O'C 
within  33°C  of  the  predictions10.  Extensive  data  bases  covering  binary 
metallic  systems  are  being  developed  in  the  U.S.  ,  France  ,A  ,  Germany  , 
England14’15,  Sweden16,  and  Japan17,18  which  can  be  employed  to  compute  ternary 
superalloy  systems.  Figures  3  and  4  present  recent  ternary  calculations  for 
the  Co-Cr-Zr  system  which  are  part  of  a  systematic  and  integrated  program, 
of  work  being  carried  out  at  the  National  Physical  Laboratory,  aimed  at 
identifying  and  characterizing  eutectic  behaviour  in  ternary  and  quaternary 
Co  +  Cr  or  Ni  +  Cr  based  alloys.  Calculations  of  the  analogue  systesxs  Co-Cr-Nb 
and  Co-Cr-Ta,  of  the  corresponding  Ni-Cr-X  systems,  and  of  some  quaternary 
modifications  to  these  systems  have  also  been  made.  The  phase  diagram  calcu¬ 
lations  have  been  made  using  the  well  established  CALPHAD  approach  using  a 
newly  developed  NPL  system  "ALLOYDATA"  which  stores  critically  assessed 
thermodynamic  and  structural  data  for  elements,  binary  and  ternary  alloys, 
retrieves  relevant  data  for  a  chosen  multicomponent  system  and  calculates  a 
wide  range  of  types  of  two  and  three-phase  equilibria  in  ternary  and  quaternary 
systems.  The  unusual  feature  of  the  sections  shown  in  Figures  3  and  4  is  the 


stability  of  the  (COjCr^Zr, 
Laves  phase  clear  across  the 
ternary  section.  The 
computed  isothermal  sections 
for  1300-1900K  indicate  the 
composition  and  temperature 
•ranges  over  which  a  number  of 
eutectic  troughs  exist,  and  also 
summarize  the  composition 
ranges  over  which  the  em¬ 
brittling  sigma  phase, 

Co  ,Cr  ,,  is  calculated  to 

.4  .0  .g 

exist  .  Of  parti¬ 
cular  relevance  to  the  aim  of 
the  calculations  is  the 
Figure  3. Calculated  1700°K  Co-Zr-Cr  Section  15  eutectic  valley  between  fee 

Co-based  solutions  (a)  and 

the  (Co.CrJjZr  Laves  phase  which  may  be  of  potential  use  for  turbine  blade 
applications.  The  isothermal  sections  indicate  the  solidification  path  along 
this  valley,  and  demonstrate  how  the  detailed  calculation  of  phase  diagrams 
can  assist  in  the  understanding  of  fairly  complex  solidification  processes. 

The  calculated  compositions  of  this  eutectic  valley  have  been  verified  by 
subsequent  experimental  studies  to  be  within  0.5-2  mol  %,  although  the 
calculated  eutectic  temperatures  are  somewhat  higher  than  the .experimentally 
determined  values.  Bearing  in  mind  the  lack  of  experimental  thermodynamic 


data  for  the  Co-Zr  system,  and 
the  uncertain  nature  of  the 
liquidus  and  solidus  curves  for 
the  Co-Cr  system,  the  agreement 
is  very  satisfactory15. 

2.  Sigma  Phase.  Due  to  the 
wide  range  of  applicability  of 
•nickel  and  cobalt  base  super¬ 
alloys  and  of  high  chromium 
stainless  steels  substantial 
effort  has  been  directed  toward 
computing  ternary  diagrams 
establishing  the  stability 
range  of  the  sigma 
phase  .  Figure  5  shows 


a  recent  result  by  Chart, 


Figure  4. Calculated  1600#K  Co-Zr-Cr  Section15 
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Figure  5.  Calculated  Stability  Range  of  the  Sigma  Phase  in  the  Cr-Fe-Ni-Si 

System  at  800K,  for  Alloys  Containing  up  to  0.08  Molt  Fraction  Si.j4 
Phase  Relationships  in  the  Ternary  Cr-Fe-Ni  Systems  are  Indicated  1 

14 

Put land  and  Dinsdale  for  quarternary  section  in  the  Fe-Cr-Ni-Si  system  at 

800K.  This  system  is  of  particular  interest  in  the  development  of  stainless 

steel  tubing  for  nuclear  reactors  which  is  resistant  to  void  swelling  under 
19  20 

irradiation  ’  .  A  new  procedure  for  evaluating  the  sigma  phase  precipitation 

tendency  of  superalloys  has  been  developed  by  Machlin  and  Shao  .  It  is  called 
SIGMA-SAFE.  Phase  diagrams  were  used  to  determine  the  degree  of  supersatura¬ 
tion  with  respect  to  precipitation  of  the  sigma  phase  as  a  function  of 
temperature  and  alloy  compositions.  It  was  found  that  this  supersaturation 
number,  A,  distinguishes  between  24  superalloys  according  to  their  known 
tendency  to  precipitate  or  not  to  precipitate  the  sigma  phase,  except  for  one 
superalloy.  This  capability  of  predicting  the  sigma  phase  precipitation 
tendency  of  an  arbitrary  new  superalloy  composition  is  not  shared  by  PHACOMP. 

It  was  shown  that  A  correlates  quantitatively  to  wt  pet  sigma  precipitated 

in  a  specific  alloy.  SIGMA-SAFE  also  yields  the  multicomponent  phase  diagram 

21 

applicable  to  each  superalloy  .  The  method  depends  upon  fixing  the  geometry 
of  the  two  and  three  phase  fields  bounding  the  y,y'  and  a  fields  in  multi- 
component  nickel  base  superalloys  as  illustrated  in  Figure  6.  Machlin  and 

g 

Shao  carried  this  out  with  the  help  of  computed  isothermal  ternary  sections 
hav«»  aoreement  as  illustrated  Table  1. 
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TftbU  I.  Comp* ri ion  of  SIGMA  SAFE  Prediction*  with  Obtemtiom  21 
*  Alloy*  Observe d  lo  Be 


Sigma  front 

Sigma  Free 

INC0713C 

(0  243) 

INCONEL  700 

IN  100 

(Borderline) 

HARM  200 

PDRL  103 

<0!/0l) 

MAR  M  246 

IN  728  NX 

(1333) 

TRW  1900 

TRW  1800 

(0.46$) 

N1CROTUNG 

U  500 

(0  393) 

TRW  NASA  IVY 

U  700 

(06S2) 

TRW  NASA  VIA 

Nimonic  115 

(Borderline) 

WASPALOY 

US20 

(0.93) 

UMTEMPAF21D 

RENE 41 

(1.3*2) 

UMTEMP  AF21DA 
NIMONIC  SO 
NIMONIC  90 
■1900 

H232 

m  Value  In  fcrcnlhcUt  I*  t  lie  supcrtaluraliun  number,  defined  to  equal  tlwe  mint* 

Flg.6— Perspective  view  of  high  Nl  region  of  Ni-Cr-Mo-Al  mum  value  of((pct  Al)  +  (pet  "((pvt  Al)  ►  (pel  Mo))o  where  the 

quaternary.  Cross-hatched  region  defines  y  ♦  y'  +  o  3  phase  arrow*  denote  the  vector*  parallel  to  ihc  (pa  Al)  and  (pet  Mo)  coordinate*  in  the 
region.  Below  and  above  tho  cross-hatched  region  the  re-  conitant  Cr  section  of  the  quaternary  corresponding  to  the  Cr  content  of  the 
•pcctlvc  3  phase  regions  each  contain  y  ♦  y'  +  a. 21  supcralky. 

.'The  assessment  of  their  SIGMA-SAFE  proposals  shows  that  for  the  10  superalloys 
known  to  be  prone  to  precipitation  of  the  sigma  phase  at  875<>CJ  SIGMA-SAFE 
correctly  predicts  all  to  be  supersaturated  with  respect  to  precipitation  of 
the  sigma  phase.  Out  of  a  total  of  14  superalloys  reported  to  have  been  held 
at  875°C  for  long  time  and  in  which  no  sigma  phase  was  found  13  were  predicted 
by  SIGMA-SAFE  to  be  in  equilibrium  phase  fields  that  do  not  contain  sigma  while 
one  was  predicted  to  lie  in  a  y+a.  field  (M252) .  The  supersaturation  number 
predicted  by  SIGMA-SAFE  correlates  to  wt  pet  sigma  quantitatively,  as  well  as 
does  the  Ny  number  predicted  by  PHACOMP.  The  critical  composition  for  preci¬ 
pitation  of  sigma  in  IN  100  is  correctly  predicted  by  SIGMA-SAFE.  The  bound¬ 
aries  between  the  y,y+a,y+y'  and  y+y’+o  phase  fields  in  the  multicompnent 
phase  diagrams  corresponding  to  Ni  base  superalloy  compositions  are  predicted 
by  SIGMA-SAFE  as  a  function  of  temperature.  The  y’  solvus  temperature  predict¬ 
ed  by  SIGMA-SAFE  agrees  with  those  observed  for  4  superalloys  within  a  100°C 
maximum  deviation.  The  predictions  of  PHACOMP  and  SIGMA-SAFE  for  the  effects 
of  specific  alloying  elements  in  promoting  sigma  differ  quantitatively,  and 
for  some  elements,  qualitatively.  This  difference  can  be  important  for 
certain  superalloys. 

3.  High  Temperature  Battery  Electrodes.  Solid  Li-Al  alloys  appear  to  be  one 

-of  the  leading  candidates  for  the  negative  electrode  of  high  temperature 

batteries.  Nevertheless  a  loss  of  capacity,  accompanied  by  changes  in 

morphology  of  the  Li-Al  components  has  been  encountered  after  cyclic  operation. 

Modification  of  cell  design  and  alteration  of  the  electrode  composition  is 

being  considered  in  an  effort  to  overcome  the  above  noted  shortcomings  and 

22 

improve  Li  diffusion  in  the  a-Al  solution  .  Addition  of  Mg  has  been 
considered  as  a  viable  solution,  however  such  a  route  requires  knowledge  of 
the  Al-Li-Mg  phase  diagram.  To  fill  the  gap  in  experimental  data,  computations 
of  the  component  binary  systems  and  the  ternary  sections  between  648  and  773K 
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Figure  7.  Calculated  Isothermal  Sections  in  the  Ca-Mg-Li  System  . 
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were  performed  .  These  computations  define  the  field  of  temperatures- 

composition  stability  for  alloys  in  this  system.  Recently  these  computations 

have  been  extended  to  cover  the  binary  components  and  ternary  sections  in  the 
.  23 

Ca-Mg-Li  system  ,  Figure  7  shows  the  results  obtained  by  Hsu  and  Saboungi 
between  673  and  773K  in  order  to  define  the  composition  range  of  stability  of 
the  solid  and  liquid  phases.  This  system  is  of  interest  in  the  development- 
of  high-temperature  batteries  in  which  iron  sulfide  is  the  positive  electrode, 
calcium-magnesium  alloys  are  the  negative  electrode  and  mixtures  of  LiCl-KCl- 
.CaClj  salts  serve  as  electrolytes.  Thus, reaction  of  the  lithium  from  the 
electrode  with  the  Ca-Mg  alloy  can  lead  to  formation  of  ternary  Ca-Mg-Al 
alloys  whose  field  of  stability  can  be  specified  by  Figure  7  for  67?-773°K. 

4.  Hydrogen  Storage  Compounds.  The  intermetal lie  compound  FeTi  is  most  likely 
to  be  the  first  to  undergo  commercial  development  for  energy  storage  purposes. 
It  can  be  hydrided  at  room  temperature  under  a  few  atmospheres  pressure;  the 
hydrogen  can  subsequently  be  released  at  a  few  tens  of  atmospheres  by  heating 
to  50-100#C.  This  range  of  pressure  and  temperature  make  iron-titanium  h; iride 
a  natural  energy  storage  medium  for  peak-shaving  systems,  solar  energy 
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Figure  8.  Recalculated  (Left  Panel)  Composite  (c)  and  Experimental  Equilibrium (d) 
in  the  Manganese-Titanium- Iron  System  at  1273K24. 


Figure  9.  Recalculated  (Left  Panel),  Composite  (c)  and  Experimental (d) 

Equilibrium  in  the  Aluminum-Titanium- Iron  System  at  1273K24.  • 

conversion  and  storage,  and  fuel  for  transportation  systems. 

The  hydriding  characteristics  of  the  FeTi  intermetallic  can  be 
manipulated  by  addition  of  a  third  element.  Replacement  of  iron  by  manganese 
improves  the  stability  of  the  compound  hydride,  making  it  more  suitable  for 
peak-shaving  applications.  It  is  important  to  have  single  phase  FeTi.  Excess 
iron  yields  Fe2Ti  Laves  phase,  which  does  not  hydride.  Excess  titanium  leads 
to  solid  solution  phase,  whose  hydride  if  so  stable  that  the  hydrogen  is  not 
recoverable.  The  presence  of  either  reduces  the  hydriding  capacity  and 
efficiency,  of  the  alloy.  Knowledge  of  the  iron-titanium-manganese  ternary 
phase  diagram  is  necessary  for  producing  single  phase  material.  In  addition, 
the  alloy  cost  is  high,  directly  melted  from  sponge  titanium  and  scrap  iron, 
the  cost  is  estimated  to  be  from. $2. 50-3. 00  per  lb.  The  compound  could  be 
produced  by  direct  reduction  from  ilmenite  (an  ore  which  contains  Fe  and  Ti 


Partltl  vertical  (action  at  0  05%Ti  for  tha  calculated  Partial  vartical  (action  at  0>50XSi  for  the  calculated 
Af-Ti-Si  ty(tam  (L— liquid.  0-TiAI,|.  Al-Ti-Si  aycam  (L — liquid.  0-TiAI,). 
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Figure  10.  Calculated  Vertical  Sections  in  the  Al-Ti-Si  System  . 

in  about  equal  proportions)  at  a  lesser  cost.  Presently  only  aluminothermic 

reduction  is  thought  to  offer  cost  reductions  below  $2.00  per.lb.  However, 

this  process  will  lead  to  aluminum  in  the  alloy.  Again  single  phase  material 

is  desirable,  and  a  knowledge  of  the  iron-titanium-aluminum  phase  diagram  is 

essential.  To  provide  such  information,  synthesis  of  the  manganese-titanium- 

iron  and  aluminum-titanium-iren  systems  were  carried  out  by  using  the  analytical 

descriptions  developed  for  the  Fe-Ti,  Fe-Mn,  Ti-Mn,Al-Ti  and  FeAl  binary 

systems  and  combining  these  to  generate  ternary  sections  \  To  provide 

additional  experimental  information  a  study  of  selected  portions  of  the  Mn-Ti-Fe 

and  Al-Ti-Fe  system  were  carried  out  at  1000°C(1273K) .  These  sections  are 

shown  in  Figures  8  and  9  with  a  recalculation  of  the  Mn-Ti-Fe  and  Al-Ti-Fe 

24 

systems  to  provide  closer  agreement  with  the  experimental  results  . 

5.  Crain  Refinement  of  Aluminum.  Grain  refinement  of  A1  by  Ti  increases 

markedly  when  the  Ti  concentration  exceeds  0.08  a/o  resulting  from  the  onset 

of  a  peritectic  reaction.  Evaluation  of  the  role  of  Si  additives  has  been 

25 

carried  out  by  Youdelis  by  constructing  the  vertical  sections  shown  in 
Figure  10  from  computed  Al-Ti-Si  isothermal  sections  generated  between  700 
and  1000K  from  binary  descriptions.  These  results  show  that  silicon  additions 
tend  to  move  the  liquidus  line  defining  the  two  phase  (L+6)  field  to  lower  Ti 
concentrations  thus  lowering  the  critical  Ti  concentration  (for  l'ie  onset  of 
the  peritectic)  to  lower  levels.  This  calculated  result  agrees  with  the 
observation  that  enhanced  grain  refinement  is  observed  at  silicon  levels  as 
low  as  0.05  a/o. 
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figure  11.  Different  types  of  three-phase  equilibrium  between  ct, 3  and  L  phases 
in  relationship  to  the  position  of  the  intersection  point  of  the 
(xo,To)  curves26. 

6.  Metatectic  Reactions  in  Hafnium  Alloys.  The  netatectic  reaction  in  which 

(Ho+L  on  cooling  has  a  deleterious  effect  on  the  high  temperature  properties 

of  alloys.  This  transformation  can  cause  excessive  segregation  during 

solidification,  welding  or  heat  treatment.-  Moreover,  when  there  is  limited 

solubility  between  low  temperature  phases,  a  eutectic  or  peritectic  results 

at  lower  temperatures  ,  which  reduces  the  capability  of  the  alloy  for  high 

26 

temperature  service.  Pieraggi  and  Dabosi  have  used  a  regular  solution  model 
to  predict  the  occurence  of  this  reaction  in  hafnium  base  alloys  as  illustrated 
in  Figure  11.  In  particular  they  have  related  the  tendency  for  metatectic 
formation  to  the  difference  between  melting  and  transformation  temperatures  of 
the  solvent  element  and  to  the  strain  energy  contribution  to  the  interaction 
parameters. 
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7.  Hardenability  of  Alloy  Steels.  Kirkaldy,  Thomson  and  Baganis 


coupled  thermochemical  and  kinetic  data  for  alloy  steels  to  develop  an 
efficient  method  for  computing  CCT  and  TTT  diagrams  illustrated  in  Figures 
12  and  13.  This  has  been  carried  out  by  using  free  energy  changes  and  activity 
data  for  iron  and  its  binary  and  ternary  alloys  which  are  used  to  evaluate  the 
general  linear  series  (Wagner)  expansion  of  the  activity  coefficient.  The  latter 
have  been  employed  in  turn  for  the  accurate  thermodynamic  determination  of  the 
liquidus  line  and  the  y  field  in  steels  with  additions  of  Mn,  Si,  Ni,  Cr,  Mo, 
Cu,  V,  Kb,  W  and  Co.  A  computer  program  based  on  several  derived  analytic 
formulae  is  accurate  for  total  additions  up  to  5  ut.%.  Indeed,  the  predictions 
compare  favorably  with  observations  on  over  200  steels  from  international 
compendia.  The  predicted  Ae^  temperatures  and  diffusion  data  have  been 
incorporated  into  a  Zener-type  expression  for  the  transformation  kinetics  of 
ferrite  and  a  semi-empirical  formula  has  been  optimized  to  adequately  predict 
the  ferrite  start  line  of  the  TTT  diagrams  of  low  alloy  steels.  The  TTT 
diagrams  are  converted  in  turn  to  CCT  diagrams  via  the  Avrami  rule.  Compari¬ 
sons  of  the  predictions  with  the  CCT  curves  for  a  representative  set  of 

27 

commercial  steels  show  good  agreement  . 

8.  Growth  of  HEM  Silicon  Solar  Cell  Material.  The  heat  exchanger  method  (HEM) 


developed  to  grow  large  sapphire  crystals  (25  cm  diam.  x  25  cm.  high)  has  been 
applied  to  the  growth  of  silicon  solar  cell  materials.  In  this  method  the 
seed  is  placed  at  the  bottom  of  the  crucible  and  the  temperature  in  the  melt 
increases  upwards.  This  suppresses  convection  that  causes  temperature  and 
concentration  fluctuations  at  the  solid-liquid  interface. 

Early  experiments  indicated  that  SIC  particles  were  found  in  crystals 
solidified  by  the  HEM.  However,  even  with  the  presence  of  SiC  particles, 
large  grains  have  been  grown  with  limited  interface  breakdown  during  solidifi¬ 
cation.  This  observation  differs  from  Czochralski  (CZ)  growth  where  interface 
breakdown  due  to  SiC  is  followed  by  twin/polycrystalline  growth. 

The  basic  elements  of  the  HEM  and  the  CZ  growth  furnaces  and  the 

processes  are  quite  similar  (heaters,  crucibles,  insulation, etc. ) .  A  silica 

crucible  loaded  with  the  charge  and  set  in  a  graphite  retainer  is  placed  in 

the  furnace.  The  chamber  is  evacuated,  and  after  melting  the  charge  crystal 

growth  is  achieved.  In  the  HEM  process,  the  chamber  is  typically  evacuated 

during  growth  to  0.1  Torr.  For  CZ  growth  an  argon  blanket  is  used  and  the 

chamber  pressure  can  vary  from  10  Torr  to  1  atm.  In  early  experiments  with 

HEM  growth  it  was  quite  surprising  that  high  carbon  concentrations  were  found 

in  the  silicon  ingots.  A  thtrmochemical  analysis  was  conducted  to  assess  the 
28 

source  of  the  problem  .  This  analysis,  illustrated  in  Figure  14,  was  subsequent 
ly  confirmed  experimentally.  Analysis  of  the  equilibrium  between  C,  CO,  CO2, 
Si,  SiO  and  SiO„  over  the  temperature  and  pressure  range  of  interest  shows 
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Figure  14.  Calculated  Pressure-Tcraperoture  Fields  of  Stability  for  SiC 

Stability  During  HEM  Growth  of  Silicon  Solar  Cell  Materials  in 
the  Presence  of  Graphite28. 
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that  CO  is  present  in  the 
working  range.  The  reaction 
of  carbon  monoxide  with  sili¬ 
con  to  yield  silicon  carbide 
and  silicon  monoxide  is  the 
source  of  the  silicon  carbide. 
One  source  of  CO  is  the  reac¬ 
tion  of  Si02  crucible  with 
the  graphite  retainer.  Re¬ 
duction  of  the  CO,  effected 
by  switching  to  a  Mo  retainer 

reduced  the  carbon  concentra- 
28 

tion  in  the  silicon  . 


Figure  15.  Calculated  Pseudobina  y  system 
(GaSb)lx(AlSb)x30 


9.  Multicomponent  Semiconduc¬ 
tors.  Simple  solution  models. 


generally  based  on  regular  solutions  are  applied  successfully  to  compute  the 

phase  diagrams  for  alloy  semiconductors,  III-V  compounds  and  a  variety  of 

multi-component  systems  which 

are  important  for  electro-optical 

and  light  emitting  diode 

applications1.  The  location  of 

the  liquidus  phase  boundaries  is 

critical  to  successful  growth 

of  such  compositions  from  l bo 
29 

melt  .  Such  calculation-  also 
provide  insight  into  the  distri¬ 
bution  of  impurities  likely  to 
be  developed  in  a  crystal  grown 
from  the  melt.  Recently, 

Linnebac^has  considered  the 
effect  of  using  a  partially 
associated  solution  (PAS)  model 
in  defining  the  solutions  in  the 
pseudobinary  system  GaSb-AlSb. 

The  results  of  these  computa¬ 
tions  are  shown  in  Figure  15. 

10.  Sulfide,  Molten  Salt,  Oxide 
and  Ceramic  Systems.  During 
the  past  ten  years,  substantial 
progress  has  been  made  in 
Jx-vcloping  data  bases,  models 


Figure  16.  The  calculated  Fe-S  diagram  with 
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isobar  values  of  S^Cg)  up  to  1  atm  . 


HQ!  and  computational  methods 

for  calculating  the  phase 
diagrams  and  thermochemical 
properties  of  metallic  sys¬ 
tems.  The  development  of 
similar  data  bases  and 
models  for  dealing  with  non 
metallic  systems  has 
proceeded  at  a  .lower  pace. 
Nevertheless,  measurable 
advances  have  been  made  and 
it  is  likely  that  the  next 
few  years  will  see  even 
.  SrCI,  more  rapid  development  of 

Figure  17.  Calculated  Liquidus  Contours  in  the  non-metallic  field 

KCl-SrCl^-CaClj  System3* .  than  the  ’’metallic"  field. 

This  is  due  to  the  fact 

that  the  "non-metallic"  group  covers  a  very  broad  spectrum  of  materials  where 
development  of  accessible  data  bases  and  models  for  estimating  properties  can 
be  quickly  translated  into  new  solutions  for  important  industrial  problems. 
Indeed  the  development  of  data  bases  for  the  "non-mctallic"  systems  has  begun 
and  prior  CALPHAD  symposia*'5  as  well  as  current  publications3140  provide 
a  measure  of  this  development.  Figure  16  shows  a  recent  example  of  a  computed 
Fe-S  phase  diagram  developed  by  Sharma  and  Chang  .  Chang  and  Hu  have  been 
developing  a  description  of  Metal-Metal-Oxygen  and  Metal-Metal-Sulphur 
systems  covering  associated 
solution,  defect  thermodynamic 
and’Wagncrl'models  for  dealing 
with  these  "nixed"  metal-non- 
metal  systems.  Similar  studies 
have  been  carried  out  by 
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Schwerdtfeger  and  kitsch  . 

Pelton,  Bale  and 
c. -ers  3^have  been 
.cooping  data  bases  cover¬ 
ing  molten  salt  and  oxide 
systems.  Figures  17  and  IS 
show  some  of  their  recent 

results35  for  fused  salt  Na<_i  KC! 

systems.  Th_  thcrnod>nnmic  Figure  18.  Calculated  l.iquidus  Contours  in  the 
properties  of  all  phases  in  LiCl-KCl-NaCl  System3*, 

the  fifteen  binary  systems 


Figure  18.  Calculated  l.iquidus  Contours  in  the 
LiCl-KCl-NaCl  System3*. 


can  be  expressed  as  differences 
of  single-ion  site  preference 
energies.  The  site  preference 
energies  and  redox  energies, 
used  to  give  the  best  agreement 
with  the  high  temperature  phase 
diagrams  are  in  reasonable 
agreement  with  values  predicted 
by  earlier  analyses  of  thermo¬ 
dynamic  and  spectral  data  as 
well  as  with  calculations  from 
crystal  field  theory.  In  the 
case  of  the  Fe^O^-Co^O^  spinel 
phase,  it  was  necessary  to 
propose  a  small  additional 
positive  regular  solution  term 
to  obtain  quantitative  agree¬ 
ment  between  the  measured  and 
calculated  phase  diagrams  at 
all  compositions.  In  the  other 
systems,  the  simple  model  with 
no  additional  terms  was  found 
to  be  sufficient36. 

Additional  efforts38'43 
have  been  directed  toward 
development  of  data  bases 
capable  of  computing  ternary 
and  higher  order  phase  diagrams 
from  the  thermochemical  proper¬ 
ties  of.  oxide,  nitride  and 

mixed  systems  which  have  direct  6  "0.5  0.5  0.4  0.6 

38 

application  in  slag-metal  AIq  ^0^  ^  System  . 

equilibria,  high  temperature 

refractories  and  (SiC^-A^Oj-SijN^-AlN)  synthesis  of  SIAL0N  ceramics.  Figure 
20  shows  an  example  of  calculated  and  observed  isothermal  sections  in  the  MgO- 
CrjOj-A^Oj  system.  .One  of  the  chief  advantages  which  this  effort  will  provide 
is  the  ability  to  generate  metastable  equilibria  data  which  can  provide  infor¬ 
mation  on  supercooled  liquids  and  the  nature  of  the  phase  diagram  when  one 
or  more  of  the  competing  phases  is  suppressed.  Such  information  is  of  groat 
practical  value  when  dealing  with  oxide  ceramics  and  glass  forming  systems. 
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Figure  19.  Calculated  Phase  Diagram  for  the  Fe-Co-0  System  in  Air 
are  Experimental. 
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formed  among  the  six  salts  LiCl,  NaCI,  KC1,  MgCl2,  CaCl2»  and  SrCl2  were 
obtained  and  expressed  in  mathematical  form  by  a  critical  analysis  of  the 
measured  binary  phase  diagrams  and  of  published  binary  thermodynamic  data. 
Interpolation  precedures  were  then  used  to  estimate  the  thermodynamic  proper¬ 
ties  of  the  twenty  ternary  systems  formed  among  the  six  salts,  and  all  twenty 
ternary  phase  diagrams  were  calculated  using  a  computer  program  designed  for 
the  purpose.  Good  agreement  with  the  measured  phase  diagrams  for  the  eight 
ternary  systems  for  which  measured  diagrams  are  available  was  obtained  3 . 

Synthesis  of  thermochemical  properties  and  phase  diagram  data  has 
also  been  performed36,37  for  the  Fe-Cr-0,  Fe-Ni-0,  Cr-Ni-0,  Mn^O^-COjO^, 
FejO^-Mn^  and  Fe^O^-Co^O^  systems.  Figure  19  shows  an  example  of  a  section 
in  the  latter  system  computed  at  an  oxygen  pressure  corresponding  to  air  at 
one  atmosphere.  In  the  former  study  high  temperature  phase  equilibrium  data 
was  compiled  and  assessed  to  evaluate  the  thermodynamics  and  structural  models 
of  the  spinel  phases.  The  thermodynamics  of  the  spinel  phases  was  determined 
in  terms  of  the  distribution  of  the  cations  between  the  tetrahedral  and  octa¬ 
hedral  sites.  The  phase  diagram  analysis  approach  to  examining  this  cation 
distribution  has  the  advantage  that  a  true  high  temperature  equilibrium  is  being 
examined.  The  high  temperature  phase  diagram  data  was  shown  to  be  consistent 
with  the  assumption  that  the  equilibrium  constants  fox  the  exchange  and  redox 
reactions  are  independent  of  composition  and  that  the  free  energies  of  the 
three  site  exchange  reactions,  but  not  of  the  redox  reaction,  are  independent 
of  temperature.  Furthermore,  the  free  energies  of  the  site  exchange  reactions 
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Conclusion 

The  foregoing  set  of  examples  show  how  far  the  utilization  of  data 
bases  and  computer  methods  for  solving  industrial  problems  has  progressed  in 
just  the  past  five  years.  Not  too  long  ago  it  was  a  difficult  task  to  consider 
computing  a  ternary  system.  Moreover,  there  were  very  few  systems  for  which 
sufficient  data  existed  which  would  permit  a  realistic  computation  to  be 
carried  out.  By  contrast,  the  current  situation  is  such  that  computation 
techniqu<  s  for  routine  calculation  of  ternary  systems  (at  reasonable  cost)  are 
available  and  are  being  applied  to  a  very  wide  spectrum  of  real  problems. 

Higher  order  systems  are  also  being  computed  where  necessary.  Moreover,  data 
bases'  covering  a  wide  range  of  materials  systems  are  developing.  The  community 
of  effort  which  has  achieved  this  progress  must  be  continued  and  attention 
should  be  paid  to  educating  the  technical  and  scientific  community  at  large 
to  the  efficacy  of  this  development. 
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VI.  Thermochemistry  of  Doped  Germanium  Borosilicate  Glasses 
for  Communications  Fibers 

Recent  efforts  at  RADC/ET-Hanscom  Field,  Massachusetts 
have  led  to  the  development  of  the  only  fiber  drawing  tower  in 
the  DoD  capable  of  producing  multikilometer  lengths  of  jacketed 
fiber  by  a  process  which  involves  drawing  the  150  micron  diameter 
fibers  from  a  glass  preform  which  has  been  passed  through  an 
annular  hot  zone.  The  preforms  from  which  fibers  are  drawn 
are  prepared  by  the  CVD  process  such  that  the  layers  (and  the 
distribution  of  elements  providing  a  refractive  index  variation) 
can  be  deposited  in  a  controlled  manner.  This  insures  light 
guiding  of  the  communications  signal.  One  current  objective  is 
improvement  of  the  induced  transmission  loss  due  to  exposure 
to  ionizing  irradiation.  At  present  the  evidence  in  hand 
indicates  that  doping  the  preforms  with  multivalent  ions  and/or 
ions  providing  deep  energy  band  levels  in  the  silicate  glass 
may  provide  the  best  route  to  radiation  resistance  improvement. 
The  ratio  of  oxidation/reduction  levels  of  the  various  dopant 
ions  such  as  Ce+3/Ce+^  or  Sb+3/Sb+^  as  a  function  of  oxidation/ 
reducing  atmosphere  has  proven  to  be  critical  in  some  of  the 
glasses.  The  preforms  are  fabricated  in  a  thin-walled  fused 
silica  reaction  tube  which  will  serve  as  a  cladding  material 
for  the  fiber  which  is  supported  on  both  ends  in  a  glassworking 
lathe.  While  rotating,  it  is  heated  with  an  oxygen-hydrogen 
burner  which  traverses  along  its  length.  Electronic  grade 
chemicals  (SiCl^,  GeCl^,  POClj)  are  kept  in  separate  containers 
at  controlled  temperatures  and  ultra-pure  oxygen  is  bubbled 
through  these  liquids.  The  chemical  vapors  entering  the 
reaction  tube  are  controlled  using  a  mass  flow  control  system 
incorporating  a  feedback  loop  to  control  the  desired  quantity 
of  oxygen  and  chloride  vapors.  Ultra-pure  BCl-j  vapors  from 
a  pressurized  tank  are  also  monitored  and  controlled  by  the 
system,  since  can  be  used  in  chemical  vapor  deposition  as  an 


intermediate  cladding  or  core  component.  The  mixture  of  oxygen 
and  reactants  is  passed  through  the  fused  quartz  tube  continuously, 
the  heat  for  the  reaction  being  provided  by  the  traversing 
burner.  The  temperature  of  the  tube  (1600-1650°C)  is  continuously 
monitored  and  controlled  using  an  optical  pyrometer.  In  the 
reaction  zone,  both  the  heterogeneous  reaction  at  the  tube 
wall  and  a  homogeneous  gas  phase  reaction  take  place.  The 
homogeneous  reaction  produces  fine  solid  patticles  of  oxide  glass 
product,  some  of  which  are  driven  to  the  tube  wall.  At  the 
high  temperature  employed,  they  are  fused  to  a  bubble-free 
glassy  film  as  the  burner  passes  along  the  tube.  After  the 
desired  number  of  passes  has  been  made,  tne  burner  temperature 
is  increased  in  excess  of  1800°C,  and  the  tube  is  collapsed 
to  form  a  waveguide  preform. 

Since  the  CVD  process  is  largely  empirical,  ManLabs, 

Inc.  has  initiated  application  of  the  computer  based  thermochemical 
and  data  bank  analysis  of  the  CVD  problem  in  support  of  the 
RADC  effort.  The  most  expeditious  route  to  provide  this  support 
is  to  utilize  the  methods  described  in  Sections  II-V  with 
specific  consideration  of  the  phase  diagrams  composed  of 
I^Oj-SiC^-GeO  and  additions  of  P2O5,  ^2^5  anc*  CeC^. 

The  current  effort  will  attempt  to  generate  state 
diagrams  for  the  glass  CVD  process  which  can  be  applied  to 
a  variety  of  materials  of  current  interest  in  optical  wave  guides. 
These  state  diagrams  should  specify  the  temperature-pressure- 
concentration  ranges  where  specific  products  can  be  formed 
starting  with  known  reactants.  These  systems  should  be  represent¬ 
ative  of  those  of  current  interest  in  the  doped  silicate  glass 
preform  technology.  Thus,  the  best  state  diagram  information 
may  consist  of  free  energy- temperature-oxygen  (or  chlorine) 
pressure  information  which  is  capable  of  being  combined  to  deal 
with  mixtures  or  solutions  of  components  in  the  gaseous  or  con¬ 
densed  state.  The  particular  glass  systems  of  interest  will  cover 
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the  GeC^-SiC^-l^Oj  system  with  the  following  dopants:  P2®5» 

Sb20g,  Sb202>  and  CeC^.  The  current  study  will  examine 

the  CVD  process  as  it  relates  to  fabrication  of  doped  germanium- 
borosilicate  glass  preforms  used  for  optical  communication  fibers 
with  special  emphasis  on  the  current  RADC/ET  mixing  chamber 
design  employed  to  investigate  introduction  of  a  variety  of 
dopants.  In  this  configuration  the  dopants  are  introduced  as  chloride 
compounds  which  are  decomposed  and  oxidized  to  be  finally 
deposited  as  oxides  at  1550  to  x650°C.  The  decomposition  and 
reformulation  steps  will  be  treated  separately  and  then  combined 
so  as  to  simulate  the  oxide  configuration  of  the  final  preform. 

Special  attention  will  be  given  to  the  oxidation  state  of  the 
dopants  in  the  germanium-borosilicate  glass. 

In  order  to  illustrate  how  the  current  technique  can 
be  applied  to  generate  such  info  .tion  it  is  instructive  to 
consider  the  following  examples: 

a.  Utilization  of  the  ManLabs-NPL  Data  Bank  to  Define 
Vapor  Pressure  of  Dopants 

Thermochemical  data  on  over  2000  inorganic  compounds 
are  held  in  the  data  bank.  Many  of  the  halides  of  the  dopants 
cf  interest  are  included.  The  chlorides  of  antimony,  SbClj-  and 
SbClj  in  the  gaseous  and  condensed  state, are  included.  Figure 
1  and  Table  1  illustrate  how  the  bank  can  be  employed  to  compute 
the  vapor  pressure  of  SbCl5  over  solid  SbCl^  as  a  function  of 
temperature.  Table  1  shows  the  reaction  -C15SB  +  C15SB(G),  which 
is  the  computer  language  for  the  vaporization  of  solid  SbClj. 
displayed  in  Figure  1.  Entry  of  this  notation  in  the  ManLabs-NPL 
system  causes  the  computer  to  output  the  specific  heat,  entropy, 
enthalpy  and  free  energy  change  for  the  reaction  as  a  function 
of  temperature  as  shown  in  Table  1.  In  addition  LOGIOKP  is 
tabulated.  The  latter  can  be  employed  to  compute  the  vapor 
pressure  of  SbCl^  since 

Kp  =  P(SbCl5) (atm)  =  P(SbCl5) (torr)/760  (1) 
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TABLE  1 

DATA  BANK  OUTPUT  FOR  VAPORIZATION  OF  SbCl 


73  01. 


thermcidynam  i  c  functi  arts 


•CL5SB+CL5SB  05  > 


298. 15 
300.  00 
335. 00 
350.  00 
375. 00 
400.  00 
435. 00 


BELT AGP 
CfiL/K.  MOL 


-18. 190 
-11.588 
-11. 000 
-10.572 
-10.819 
-3. 383 


DELTAS 
CflL/K  MOL 

87. 000 
86. 984 
85.976 
85, 148 
84.338 
83. 787 
83.117 


DEL T AH’ 
CAL/NQL 

11900 
1 1 37? 
11581 
11300 
11030 
10771 
10519 


DEL TAG 
CAL/riOL 

3850 
3300 
3139 
8500 
1 33 1 
1830 
694 


LOGIOKP 


“8..  388 

-8. 763 

~8. 1 1 1 
-1.56! 
-1. 096 
-0.699 


The  vapor  pressure  of  SBCl^  in  torr  calculated  from  Eq.(l) 
and  the  data  in  Table  1  is  displayed  in  Figure  1.  This  shows 
that  at  350°K(i.e.  77°C)  the  vapor  pressure  of  SbClj.  is  20.8 
torr  while  at  300CK  (i.e.  27°C)  it  is  only  1.6  torr  which  is 
too  low  for  prac  :ical  use.  In  fact,  the  20.8  torr  pressure 
is  used  by  maintaining  the  SbCl^  feed  lines  at  77°C. 

b.  Calculation  of  the  Relative  Proportion  of  SbCl^ 
and  SbClg  as  a  Function  of  Temperature  and 

Pressure 

The  foregoing  example  illustrates  how  the  data  bank 
can  perform  a  simple  problem.  A  second  example  of  a  more 
complex  problem  can  be  illustrated  by  considering  the  decomp¬ 
osition  of  SbCl5  into  SbCl^  and  Cl2  which  can  be  written  as 

SbCl3  +  Cl 2  ♦  SbClr  (2) 

Table  2  shows  the  Data  Bank  results  for  the  forward 
reaction.  The  fifth  column  tabulates  AG°  for  the  forward 
reaction, thus  at  equilibrium 

AG=0=AG°  *  RTln  /(Pc,  ^  CPsbclj)  W 

Since  PC1  _  PSbC1  =  P  x  (where  is  the  total 
2  3  t  t 

pressure  and  x  is  the  fraction  of  SbCl^)  and  =  P^Cl-Zx). 

Equation (3)  can  be  recast  as 

-AG°=RTln  (l-2x)/x2Pt 

thus  permitting  solution  of  x  (and  l-2x)  as  a  function  of 
temperature  and  pressure.  Figure  2  displays  the  results  for 
P  *1  atm  and  Pt=0.1  atm  (76  torr).  The  results  show  that 
at  1  atmosphere  and  300K  the  mixture  consists  largely  of  SbClg 
with  very  little  SbClj.  As  the  temperature  is  increased  the 
SbClj  fraction  increases  substantially.  Lowering  the  pressure 
to  76  torr  has  a  very  marked  effect  on  the  ratio.  Thus  control 
of  the  pressure  could  be  an  easy  way  to  vary  the  ratio  of 
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TEMPERATURE,  T°K 

Figure  2.  Calculated  Fraction  of  SbCl3  and  SbCl5  as  a  Function 
of  Temperature  and  Total  Pressure  in  Chlorination  of 
SbCl3 


TABLE  2 

DATA  BANK  OUTPUT  FOR  CHLORINATION  OF  SbCl3 


:L2  03)  — S 

BCL3  u3>  +SJBCL5 

i.  3) 

T 

DELTflCP  • 

DELTR2 

DELTRH 

DELTRb 

L  3b 1 0KP 

K 

CriL/K  MOL 

CRL/K  MOL 

CRL^MOL 

CflL/MDL 

293. 15 

0.  170 

—35. 290 

-13450 

—7928 

5.312 

300.  00 

0.1 99 

-35. 289 

-13450 

-7363 

5.728 

325. 00 

0.550 

-35. 259 

-18440 

-6981 

4.695 

350.  00 

0.330 

-35.207 

-18423 

-6100 

3.809 

375. 00 

'  1.057 

—35. 142 

-18399 

-5221 

3.  043 

400.  00 

.  1.244  ' 

—85. 063 

-18370 

-4343 

2.373 

425: 00 

1.400 

—34-987 

-13337 

—3463 

1 . 783 

450.  00 

•  1 . 532 

-34.903 

-18300 

-2594 

1.260 

475. 00 

1.645 

—34. 818 

-13261 

-1722 

0.792 

500.  00 

1.742 

-34.731 

-13213 

-353 

0.373 

525. 00 

1 . 827 

-34.644 

-13174 

14 

-0.  006 

550.  00 

1.901 

-34.557 

-13127 

379 

-0.349 

575. 00 

1.967 

—34. 471 

-18079 

1742 

-0.  662 

600.  00 

2.  026 

-34.336 

-18029 

•  £603 

-0.948 
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+  3  +5 

Sb  /Sb  ions  entering  the  reaction  chamber. 

Thus  the  data  bank  offers  a  means  for  investigating 
+  3  +5 

the  stability  of  Sb  /Sb  interchanges  in  the  reaction 

6SbCl5(G)  +  5Sb203(S)  -  10SbCl3(G)  +  3Sb205(S)  (5) 

as  a  function  of  temperature  and  pressure.  Similar  investigation 
can  be  undertaken  for  chloride  and  oxides  of  phosphorous,  and 
cerium  in  order  to  evaluate  their  tendency  to  yield  a  pre¬ 
dominate  valence  state  in  the  glass. 

Although  the  CVD  process  is  controlled  by  kinetic 
factors  as  well  as  by  thermodynamic  factors  and  the  former 
may  prevent  attainment  of  equilibrium,  the  equilibrium  calcu¬ 
lations  can  still  serve  as  a  useful  guide  in  outlining  the 
field  in  which  certain  species  predominate  (i.e.  Figure  2). 

Since  the  dopants  finally  are  deposited  in  the  glass 
phase  it  is  proposed  that  the  techniques  developed  in  the 
present  program  for  computing  multicomponent  oxide  phase  diagrams 
be  applied  in  dealing  with  the  condensed  state  reactions  which 
take  place  during  the  preform  synthesis  at  high  temperatures. 

As  the  first  step  in  this  ana'ysis  the  Ge02-Si02, 
B203-Si02  and  B203-Ge02  binary  systems  have  been  analysed  along 
the  lines  described  in  Sections  II-IV.  The  results  are  shown 
in  Table  3  and  Figures  3-5.  The  binary  solution  phases  in  these 
systems  are  described  by  equations  similar  to  Equation  1  on 
page  27.  These  results  can  be  employed  to  compute  ternary 
sections  of  the  B203-Ge02-Si02  system  between  1600K  and  1000K 
which  are  shown  in  Figures  6  and  7.  The  latter  are  computed 
by  employing  the  des  .ription  contained  in  Table  3  and  ternary 
solution  phase  free  energies  similar  to  those  depicted  by 
Equation  4  on  page  40  with  a  zero  valued  ternary  interaction 
parameter  for  each  phase.  The  current  program  will  next  consider 
effects  on  this  system  due  to  additions  of  P20,-,  an<^  ^2^5* 
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TABLE  3 


SUMMARY  OF  LATTICE  STABILITY  AND  SOLUTION  PHASE  PARAMETERS 
FOR  GERMANIUM  BOROSILICATE  SYSTEMS 

(All  Units  in  Joules  per  gram  atom  (mole  of  atoms), T  in  Kelvins) 
BO  =  (1/5) B203  ,  SO  =  (l/3)Si02  ,  GO  =  (l/3)GeOz 


structure»  R=low  temperature  Ge02  structure,  X=Crystobalite , 
T=Tridymite,  H=a  quartz,  B=B  quartz  structure,  U=low  temperature  GeO 


GOGOLH* 


(Liquid) - (1/3) Ge02  (a  quartz) 


(l/3)Ge02 

SOSOLQ=  -6.15T 
SOSOHU- -4791-5. 69T 
GOGOLH=  3891-2. 80T 
G0G0HU-  7443- 5. 69T 
GOGORU=  7109-5. 40T 
GOGOLX=  -1.67T 
GOGOLT=  -2.01T 
GOGOLQ=  -6.15T 


BOBOLQ=4448-6 . 15T 
BOBOLX=  -1.67T 
BOBOLT=  -2.01T 
BOBOLH=  -2.80T 
BOBOLR=  -3.10T 
BOBOLU=  -8.49T 


LBOSO=64852-41 . 84T,  LSOBO=- 4184  +  20 . 92T 
XBOSO=XSOBO-TBOSO=TSOBO=HBOSO=HSOBO= 83680 
RBOSO=RSOBO=QBOSO=QSOBO=33472 


LGOSO=LSOGO=HGOSO=HSOGO=RGOSO=RSOGO=4184+12 . 552T 
TGOSO=TSOGC=XGOSO=XSOGO=83680 ;  UGOSO=USOGO=41840 

LBOGO=LGOBO=12 552  ,  QBOGO=QGOBO=33472 ,  HGOBO=HBOGO=83680 
RBOGO=RGOBO=83680 


*  These  differences  specify  the  free  energy  of  one  phase 
(i.e.  liquid)  minus  the  free  energy  of  the  second  phase 
(i.e.  for  a  given  compound). 
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Figure  7.  Calculated  Isothermal  Sections  in  the  B203  -  Ge02-Si02  System 


